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SUMMARY 


The  vertical  gust  response  and  its  alleviation  for  hingeless  helicopter 
rotor  blades  in  cruising  flight  is  studied  theoretically  and  experimentally. 
An  evaluation  is  performed  of  the  effectiveness  of  torsional  stiffness  varia 
tion  in  conjunction  with  chordwise  center-of-gravity  shift  in  alleviating 
the  blade  flapping  response  to  decrease  the  root  bending  moment. 

The  theoretical  analysis  utilizes  the  equations  of  motion  of  hingeless 
rotor  blades  exposed  to  vertical  gusts  in  forward  flight  for  the  flapping, 
lagging,  and  elastic  and  rigid  pitch  degrees  of  freedom.  The  equations 
include  the  effect  of  steady-state  deflections  in  the  trim  conditions  and 
various  hingeless  rotor  configurations  such  as  precone,  droop  and  torque 
offset  as  well  as  chordwise  center-of-gravity  shift  and  aerodynamic  center 
offset. 

The  experimental  program  involves  the  wind  tunnel  tests  of  a five-foot 
diameter  rotor  subject  to  a sinusoidal  waveform  gust.  Testing  involves 
variation  of  the  blade  chordwise  center-of-gravity  location,  the  blade 
torsional  stiffness,  rotor  advance  ratio,  and  vertical  gust  frequency. 
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SECTION  1 


INTRODUCTION 


1.1  General 

Over  the  last  decade  the  hingeless  rotor  has  attained  acceptance  as  an 
attractive  concept  for  rotorcraft,  and  considerable  research,  development,  and 
testing  has  been  conducted.  The  main  advantages  of  this  system  are  improved 
flying  qualities,  and  design  simplicity  of  the  rotor  and  rotor  hub.  The  improved 
flight  characteristics  are  achieved  by  transmission  of  the  large  blade  moments 
to  the  aircraft  because  of  the  cantilevered  boundary  condition  at  the  blade 
root  of  the  hingeless  rotor,  while  the  conventional  articulated  rotor  has  re- 
duced control  power  and  damping  due  to  the  flapping  hinge  pin  which  eliminated 
the  large  bending  moment.  Higher  control  power  and  angular  motion  damping 
are  thus  obtained  in  the  hingeless  rotor  flying  qualities  and  the  maneuver- 
ability are  greatly  improved.  Another  advantage  is  the  simplicity  of  the  rotor 
and  the  hub,  which  leads  to  lower  manufacturing  cost,  greater  weight  and  space 
savings,  and  simplified  maintenance. 

However,  rotor  blade  aeroelastic  stability  and  response  become  more  sensi- 
tive to  the  various  hingeless  rotor  configurations  because  of  structural  coupling 
between  bending  (flapwise  and  chordwise)  and  torsion  of  the  long  slender  canti- 
lever blade.  Also,  hingeless-rotor  helicopters  without  feedback  systems  or  with- 
out large  horizontal  tails  have  weak  longitudinal  dynamic  stability  with  respect 
to  attitude  in  the  upper  speed  range  and  are  sensitive  to  atmospheric  turbulence. 

Therefore,  the  study  of  the  response  of  the  hingeless  rotor  to  atmospheric 
turbulence  is  necessary  in  evaluating  the  blade  fatigue  characteristics  and  the 
riding  comfort  of  the  vehicle,  and  in  the  development  of  gust  alleviation  sys- 
tems. This  study  will  aid  in  attaining  better  hingeless  rotor  flying  qualities. 
However,  like  the  fixed-wing  aircraft  for  which  extensive  literature  on  analysis 
and  testing  of  turbulence  response  exists,  very  little  information  is  available 
on  the  response  of  rotorcraft  to  turbulence. 

In  view  of  such  deficiencies  mentioned  above,  it  is  necessary  to  perform 


1 


further  theoretical  analysis  to  establish  an  analytical  method  for  determining 
hingeless  rotor  blade  response  to  atmospheric  turbulence  in  forward  flight 
and  to  develop  a gust-alleviation  system.  In  this  report,  blade  chordwise  mass 
balance  will  be  chosen  as  a principal  ingredient  to  the  gust  alleviation  sys- 
tem because  this  will  provide  a very  simple  blade-integrated  feedback  system 
equivalent  to  controlling  blade  pitch. 

1.2  Brief  Survey  of  Past  Work 

A good  review  and  physical  descriptions  of  the  various  dynamic  and  aero- 
elastic  problems  associated  with  hingeless  rotorcraft  have  been  given  by 
Hohenemser  in  Ref.  1.  Therefore,  references  pertinent  only  to  the  present 
problem  will  be  cited  below. 

It  may  be  said  that  the  linear  equations  developed  by  Houbolt  and  Brooks 
in  Ref.  2 were  the  first  comprehensive  equations  applicable  to  rotating  elas- 
tic cantilevered  beams  with  twist.  These  free-vibration  equations  have  been 
further  developed  by  Hodges  and  Dowell  in  Ref.  3 to  include  nonlinear  terms 
expressing  large  steady-state  deflections. 

In  the  early  analyses  of  hingeless  rotor  instability.  Young  (Ref.  4)  in- 
vestigated the  role  of  nonlinear  coupling  terms  in  flap-lag  instability. 
Subsequent  investigation  was  done  by  Hohenemser  and  Heaton  (Ref.  5),  using  a 
numerical  integration  scheme  to  analyze  flap-lag  instability  in  hover  and  for- 
ward flight.  A linear  stability  analysis  of  the  hingeless  blade  in  hover  was 
conducted  by  Ormiston  and  Hodges  (Ref.  6) . They  introduced  structural  flap- 
lag  coupling  to  supplement  the  spring-restrained,  centrally-hinged,  rigid  blade 
model  which  most  predecessors  used.  The  perturbation  method  in  multiple  time 
scales  was  applied  to  the  nonlinear  flap-lag  instability  problem  of  the  spring- 
restrained,  centrally-hinged,  rigid  blade  by  Tong  (Refs.  7 and  8) . The  next 
logical  step  was  to  study  the  flap-lag  instability  of  elastic  cantilever  blades. 
Friedmann  and  Tong  (Refs.  9 and  10)  studied  flap-lag  instability  in  hover  and 
forward  flight  using  the  perturbation  method  in  multiple  time  scales  and  a 
numerical  integration  method.  Friedmann  extended  his  work  in  Refs.  11  and  12. 
Peters  (Ref.  13)  studied  the  flap-lag  instability  in  forward  flight  and  showed 
the  validity  of  the  linearized  perturbation  equations  for  y < 0.5. 
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Concurrent  with  these  flap- lag  instability  analyses,  flap-lag-torsion 
equations  of  the  hingeless  rotor  blade  were  developed  and  investigated. 

Miller  and  Ellis  (Ref.  14)  developed  the  flap-lag-pitch  equations  and  showed 
the  results  of  the  effect  of  distributed  torsion  and  finite  steady-state  de- 
flections. 

Friedmann  and  Tong  (Ref.  10)  investigated  the  elastic  flap-lag-pitch 
equations.  For  the  torsion  degree  of  freedom,  the  rigid  pitch  due  to  control 
link  flexibility  was  chosen.  Friedmann  extended  the  analysis  to  incorporate 
elastic  torsion  motion  (Ref.  15) . 

Hodge  and  Ormiston  applied  their  flap-lag-torsion  equations  to  the  uni- 
form cantilever  beam  in  hover;  they  extended  their  work  in  Ref.  16.  In  Ref.  17 
they  incorporated  rigid  pitch  motion  due  to  control  link  flexibility. 

Extensive  investigations  were  conducted  on  the  hingeless  helicopter 
BO-105  using  the  rigid  offset-hinged  blade  mathematical  model  (Refs. 18  through 
24).  They  described  various  important  design  considerations  for  the  hingeless 
rotor  and  showed  analytical  results  in  good  agreement  with  experimental 
data.  They  also  showed  the  significance  of  structural  coupling  such  as  pitch- 
lag  coupling.  Hansford  and  Simons  (Ref.  25)  also  mentioned  the  important  role 
of  bending-torsion  coupling. 

Unfortunately,  very  little  literature  is  available  on  the  response  of  the 
helicopter  to  atmospheric  turbulence  as  mentioned  before.  Among  those.  Refs. 
26,  27,  28,  and  29  analyzed  the  random  blades  response  to  random  inputs,  using 
various  methods  of  computing  the  random  response  of  a time-varying  linear  sys- 
tem. Gaonker  and  Hohenemser  (Ref.  30)  also  showed  the  adequacy  of  the  assump- 
tion that  the  turbulence  velocity  is  uniform  over  the  entire  rotor  disk  if  the 
ratio  of  the  turbulence  scale  length  to  the  rotor  radius  is  larger  than  4.0. 


Arcidiacono  et  al.  (Ref.  31)  conducted  an  analytical  study  of  articu- 
lated helicopter  response  to  various  gust  profiles  such  as  sine-squared  and 
ramp  gusts.  They  confirmed  that  the.  current  gust  alleviation  factors  speci- 
fied in  MIL-S-8698 (ARG)  are  too  conservative.  However,  they  pointed  out  that 
the  gust  loadings  can  be  an  important  consideration  from  a ride-comfort  stand- 
point. 
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For  the  tilt-rotor  aircraft,  extensive  analyses  of  the  gust  response 
in  cruising  flight  were  conducted  by  Johnson  (Ref.  32)  and  Yasue  (Ref.  33). 

In  Ref.  32,  Johnson  summarized  the  comparison  of  theory  with  the  results  of 
full-scale  tests  of  two  proprotors  in  terms  of  critical  damping  ratios  and 
frequencies  of  modes.  The  agreement  between  the  theory  and  the  experiment  was 
good.  He  extended  his  theory  (Ref.  34)  to  include  various  rotor  configurations 
such  as  precone  and  droop  as  well  as  the  transition  flight  regime  between  hover- 
ing and  cruising  flight.  Concurrent  with  the  extension  of  the  theory,  Frick 
and  Johnson  (Ref.  35)  applied  modern  control  theory  to  develop  a gust  allevia- 
tion system  for  the  tilt- rotor.  Wind  tunnel  experiments  were  conducted  at 
M.I.T.,  using  a one-ninth-scale  tilt-rotor  model  subjected  to  vertical  and  longi- 
tudinal sinusoidal  gusts  as  well  as  collective  and  cyclic  pitch  control  inputs 
(Ref.  36) . The  experimental  results  showed  good  agreement  with  the  theoretical 
analysis  of  Ref.  33.  Further  investigations  were  performed  to  develop  a gust 
alleviation  device  for  the  gimballed  tilt-rotor  (Ref.  37),  and  subsequently, 
wind  tunnel  testing  was  conducted  to  verify  the  theory  predictions  (Ref.  38) . 

Finally,  the  chordwise  blade  mass  balance  system  is  briefly  reviewed. 

Miller  (Ref.  39)  suggested  chordwise  blade  mass  balance  to  improve  the  heli- 
copter flying  quality.  Hirsh  et  al.  (Ref.  39)  applied  chordwise  mass  balance 
to  the  articulated  rotor  helicopter  XH-17  and  reduced  the  vibration  stresses 
from  35  to  45  percent.  Miller  and  Ellis  (Ref.  14)  analyzed  the  effectiveness 
of  chordwise  mass  balance  to  reduce  the  vibratory  root  shear  of  the  articulated 
helicopter.  According  to  Refs.  19  and  23,  a 2.5%  forward  shift  of  the  center- 
of-gravity  of  the  BO-105  b. ade  reduced  the  rotor  angle-of-attach  instability 
time-to-double  amplitude  from  7 to  22  seconds.  In  conjunction  with  the  improve- 
ment of  the  stability  characteristics,  a favorable  gust  response  reduction  was 
also  obtained,  as  shown  in  Refs.  23  and  24.  A 5%  center-of-gravity  forward 
shift  reduced  the  flapping  response  by  about  40%  in  trimmed  flight  at  100  knots, 
compared  with  no  chordwise  mass  balance.  Thus,  blade  center-of-gravity  forward 
of  the  aerodynamic  center  will  reduce  both  the  angle-of-attack  instability  and 
the  gust  sensitivity. 
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1.3  The  Purpose  and  Scope  of  the  Research 

The  purpose  of  this  work  is  to  investigate  the  effectiveness  of  chord- 
wise  center-of-gravity  shift  in  conjunction  of  blade  torsional  stiffness 
variation  in  alleviating  the  blade  flapping  response  of  the  hingeless  rotor 
blade  theoretically  and  experimentally. 

As  the  first  step  of  the  theoretical  analysis,  free  vibration  equations 
for  the  flap-lag-torsion-rigid  pitch  motion  of  the  rotating  nonuniform  beam 
will  be  developed  including  precone,  droop,  torque  offset,  offsets  between 
the  elastic  axis  and  cross-sectional  center  of  gravity,  built-in  twist,  and 
collective  pitch.  The  rigid-pitch  motion  is  due  to  control  linkage  flexibil- 
ity. Based  on  these  equations,  the  natural  frequencies  and  corresponding 
normal  coupled  mode  shapes  of  a rotating  nonuniform  rotor  blade  will  be 
obtained  using  the  subspace  iteration  method  of  the  finite-element  method. 

By  the  use  of  Galerkin's  Method,  the  equations  of  motion  for  the  blade  with 
aerodynamic  forces  will  be  converted  to  modal  equations.  The  equations  will 
include  the  effect  of  steady-state  deflections  and  the  quasi-steady  assumption 
will  be  employed  for  the  aerodynamic  forces.  The  coefficients  of  the  equations 
are  periodic  in  time  in  forward  flight,  so  the  harmonic  balance  approxi- 
mation method  (Ref.  40)  will  be  used  to  avoid  the  complexity  of  periodic  co- 
efficients. These  modal  equations  with  constant  coefficients  will  be  the 
basis  for  a study  of  the  hingeless  rotor  response  to  the  vertical  gust  in 
cruising  flight. 

Eigenvalue  analysis  will  be  used  to  examine  the  rotor  system  stability 
with  the  blade  mass  balance,  because  the  hingeless  rotor  system  stability  is 
very  sensitive  to  the  rotor  configuration.  To  measure  the  effectiveness  of 
blade  chordwise  mass  balance  to  reduce  blade  gust  response,  the  frequency 
response  of  the  blade  motion  to  vertical  gusts  will  be  analyzed. 

The  experimental  program  involves  the  design,  construction,  and  testing 
of  a five-foot  diameter  rotor.  The  rotor  design  is  such  that  the  torsional 
stiffness  of  the  blade  as  well  as  the  blade  chordwise  center-of-gravity  loca- 
tion can  be  varied  during  the  various  phases  of  the  test.  The  rotor  is 
exposed  to  a sinusoidal  waveform  vertical  gust  and  the  flap,  lag  and  torsion 
response  of  the  rotor  is  measured.  The  construction  and  validation  tests  are 
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performed  first  on  a three-bladed  rotor  employing  "stiff  torsion"  blade.  The 
tests  continue  on  a rotor  utilizing  two  counterweights  and  a single  blade.  In 
this  configuration,  a chordwise  center-of-gravity  shift  is  applied  to  the 
single  blade  in  conjunction  with  variable  blade  torsional  stiffness. 

The  experimental  results  are  compared  with  the  theoretical  predictions. 


SECTION  2 


THEORETICAL  ANALYSIS 

2.1  The  Equations  of  Motion  for  the  Rotating  Blade 

2.1.1  Rotor  Configuration 

The  hingeless  rotor  configuration  considered  in  this  report  is  described 
here.  A cantilever  rotor  blade  is  attached  to  a rigid  hub  which  incorporates 
built-in  precone,  the  inclination  of  the  pitch  change  bearing  with  respect  to 
the  plane  of  rotation  (Fig.  1) . Torque  offset  between  the  blade  pitch  axis 
and  the  center  of  rotation  to  minimize  the  steady-state  blade  stresses  is  in- 
cluded. Another  angular  offset  considered  here  is  droop,  the  inclination  of 
the  blade  outboard  of  the  pitch  change  bearing.  Although  built-in  sweep  is 
not  treated  here,  the  blade  coning  motion  about  the  blade  pitch  bearing  pro- 
duces sweep  at  certain  collective  pitch  settings. 

The  blade  outboard  of  the  pitch  change  bearing  can  be  rotated  about  the 
feathering  axis  by  moving  the  control  linkage.  Therefore,  the  control  linkage 
flexibility  will  allow  rigid-body  pitching  motion  of  the  blade  (called  rigid 
pitch  in  this  report).  Blade  flexibility  will  permit  flap  bending,  lead-lag 
bending,  and  elastic  torsion  motions.  The  present  analysis  treats  these  de- 
flections of  a twisted  blade  with  nonuniformities  of  mass,  bending  stiffness, 
and  torsional  rigidity  distribution.  Furthermore,  the  blade  cross-sectional 
center  of  gravity,  tension  axis,  and  aerodynamic  center  offsets  from  the 
elastic  axis  are  considered  because  these  are  important  factors  for  the  study 
of  the  mass-balanced-blade  gust  response. 

The  above  configurations  may  differ  for  certain  hingeless  rotor  designs. 
For  example,  the  Westland  WG-13  rotor  design  (Ref.  25) , positioned  the  pitch 
change  bearing  outboard  of  the  hub.  In  addition,  a flapwise  flexible  element 
inboard  and  a matched  stiffness  element  outboard  of  the  pitch  change  bearing 
are  incorporated  to  eliminate  blade  bending-torsion  coupling.  For  this  con- 
figuration, certain  modifications  are  necessary  in  the  present  analysis. 

2.1.2  Coordinate  Systems  and  Inertia  Loadings 

The  coordinate  system  x,  y,  and  z (Fig.  2)  is  rotating  with  a constant 
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rotor  angular  velocity  about  the  z axis  which  is  the  axis  of  rotation.  The 
x,  y plane  is  a hub  plane  with  the  projection  of  the  undeformed  blade  elastic 
axis  with  zero  collective  pitch  lying  on  the  y axis.  The  second  coordinate 
system  x',  y',  and  z'  is  defined  with  the  y'  axis  lying  on  the  blade  feathering 
axis  and  the  x'  axis  parallel  to  the  x axis.  The  third  coordinate  system,  x" , 
y",  and  z"  is  fixed  to  the  blade  root  and  as  such  locates  the  origin  at  the 
intersection  of  the  blade  elastic  axis  and  blade  root  chordline.  The  y"  axis 
coincides  with  the  undeformed  blade  axis,  the  x"  axis  with  the  chordline,  and 
the  z"  axis  perpendicular  to  the  x" , y"  plane. 

The  transformation  between  these  three  coordinate  systems  is  described 
as  follows:  the  x,  y,  z coordinate  system  is  translated  to  the  negative  x di- 
rection at  the  torque  offset  distance  e^  (positive  forward) . Then  this  trans- 
formed system  is  rotated  about  the  x axis  at  the  precone  angle  0^  (positive 
upward).  Thus,  the  y'  axis  is  determined,  which  coincides  with  the  feathering 
axis.  The  x'  axis  is  parallel  to  the  x axis  and  the  z'  axis  is  also  inclined 
at  the  angle  of  the  precone  from  the  z axis.  The  origin  of  the  system  is  next 
translated  along  the  new  y axis  at  the  distance  of  the  rigid  shank  r^  (a  portion 
of  the  hub).  The  system  is  rotated  about  the  new  x'  axis  at  the  droop  angle 
6D  (positive  upward) , then  the  y'  axis  lies  on  the  y"  axis.  This  intermediate 

system  is  again  rotated  about  the  y"  axis  at  the  pitch  angle  9 (positive  nose- 

m 

up) , which  includes  the  collective  pitch  0 , perturbation  control  pitch  input 

o 

0 , and  the  rigid  pitch  6 . This  scheme  makes  the  x' , y',  z'  axes  lie  on  the 

C R 

x",y",z"  axes  (Fig.  2). 

The  deformations  described  by  the  displacements  of  the  elastic  axis 
parallel  to  the  x"  and  z"  axes  are  v and  w,  respectively  (Fig.  3) . The  radial 
displacement  u parallel  to  the  y"  axis  is  included  in  this  report.  The  in- 
plane displacement  v is  positive  backward  with  respect  to  rotation.  The  out- 
of-plane  displacement  w is  positive  upward.  It  should  be  noted  that  v and  w 
are  displacements,  respectively,  parallel  and  perpendicular  to  the  blade  root 
chordline  which  is  inclined  at  the  collective  pitch  angle. 

The  x and  z axes  are  the  principal  axes  of  the  blade  cross  section  for 
P P 

the  shear  center  and  y is  the  tangent  to  the  deformed  elastic  axis.  The 
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inclindation  between  the  x" , z"  axes  and  the  x , z axes  defines  the  built-in 

P P 

angle  of  twist  0 and  elastic  torsion  4> » which  are  both  positive  nose  up 

(Fig.  3).  The  relation  between  the  x"  , y"  and  z"  axes  and  x , y , and  z axes 

P P P 

are  defined  by  the  rotation  about  x" , z"  axes  at  the  deflection  slope  angles 
v’  and  w' , and  rotated  about  the  y"  axis  at  the  angle  of  twist  and  torsion 
angle. 

Next,  the  position  vector  to  the  arbitrary  point  of  the  blade  cross 
section  is  defined.  This  point  on  the  cross  section  expressed  in  the  principal 
coordinates  is  as  follows: 


(1) 


This  vector  is  transformed  by  the  transformation  matrix  [T]  to  the  undeformed 
elastic  axis  system. 
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where 


e = Gtw  -v-  $ 


wf=  za: 

ar 


v'.iHl 

ar 


(3) 


The  built-in  angle  of  twist  is  expressed  by  0 ; <J>  is  the  elastic  torsion 

deflection;  and  v and  w are  chordwise  and  vertical  deflections,  respectively. 
The  radial  coordinate  from  the  pitch  change  bearing  along  the  undeformed 
elastic  axis  is  taken  as  r. 

The  distance  from  the  origin  of  the  x" , y" , z"  coordinate  system  to  the 
origin  of  the  principal  axes  of  the  blade  cross  section  is  defined  by 


= 


'U~ 

r + u. 

u r 


(4) 


where  u is  the  radial  displacement  given  by: 

u ' r 


+ (V)*j  dr 


(5) 


It  should  be  noted  that  the  radial  displacement  includes  only  the  geometrical 
shortening  effect  due  to  bending  deflections.  The  displacement  due  to  the 
tensile  forces  is  neglected  in  this  analysis. 

The  x" , y"  and  z"  systems  are  transformed  by  the  transformation  matrices 
[Tq]  and[T^]  which  are  related  to  the  collective  pitch  0^  and  the  droop  6D, 
respectively : 


where 
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The  acceleration  a of  the  fixed  coordinate  system  when  the  blade 
is  rotating  at  angular  velocity  ft  about  the  z axis  is: 


a 
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Although  double  underlined  terms  are  nonlinear  terms,  they  are  retained  to 
describe  steady-state  deflection  effects  on  the  blade  pitch  motion. 


The  acceleration  a^  in  the  undeformed  elastic  axis  coordinate  (x",  y" , 
z"  coordinate  system)  is  obtained  by  the  transformation. 
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- ocro  C.0 ~ 40)  - XT*  >xXr.  0 - Uj  Oro  © 

i>j~  - G^xr  - 2 Q ( pp  4-  ^DO3O0^)  xr-  Q-U^ 

*”  — ^/o4y^)  CSj>+  ©tM.)  + C — e.  + 

V So  A^O^J  xx^v.Qm  4 Xl'A>^*-0n  curoQ^ 

**?[-  (©r  4 4>)  euro  © - 2 O C~  ^ co-o  © 

■+  U_T  Ac '*’“©»<.  A^c ©3  — .O-  A^-^-Qh^  C<rO  (0^4  fljj 

^ ^P  C~  A'-^-Q  4 2.0  ( XT /A^*'-0»t.  A'-*'1-© 

— Cl  A^~- A*-^ -fGjj 


kT7 


(13) 


The  distributed  inertia  forces  P£A  due  to  the  acceleration  in  the  x",  y", 
z"  coordinate  system  are 

~ Tx 


p«=-J 

-'A 


/g  (2S4  cfc?p 


• • * * • 

hit  + m rS o 0*  + >*t  uj0p  - i>»x  & u 

+ 2-vn©Kij-  4-  ^Xl(^+^«t30juT 

-*'fl,(«+'1pfA0,  + c*>e„  (-  e„ 
v^So  a2*«-0*v  4 <u~  c/K)  ©«_  -f  kj"  a*4c.0-«.)^ 

4 'H.  - (0*4  0)  /*c*c0  4 2..Q.  c<n3  0^  ('V-'ajxjQ 


15 


; 


— k_t  ' 


D 0 - ccro  0^  oera  0 0^  + Q)  j 


rcn.  LX 


ii  4 2 yv\_  -O-  ( TJ"  cxtoGh^  + txT  A^\Qy 


4-  y Qfg  oaro  0m.  + ©*  UT'ocroQm.^ 

- kG?  [ <v0  4 r)  + Sd  a^-Qh  (-  e»  + v~Si> 

4-  C 'VT - vJ“ 00-0  G>n_)(  (^p  ^ Sd  c'°r°®  l*-Vj 
+ /»\  OCx  — LT  'c^jcQ  4 {j'/jCn.© 

-2D.  (Qk  +(£)  A^-(Q^-he)  -^{S'd^-O 

+ Sd  /C-^0,*.  ccro(0^-f  0)  4 Pp  /o^^(Q>^  + ©) 

**  LT 1 coo  4-  uj~  0 ^ 


'W.  l*T  - 2 M.H.  (.  £>p  4 Sictro©^)  ar 
- 'w-v .Cl  [ — OC^d 4-  (3pcoo0^)  +(-e0 

4 rS^>  /Vovo.0^  4 CcocQ^  4 
4 ™ Xr  [ ~(©„4  4>)  c^ro  Q - u-  0R 
“ 2 Cl  | - XT1  #^0^  CcrO  © 4 Ul/^Cn.0Ht  /dX^-0  j 

— 2 'O'  ©R  - O.  /4M-^-O0h  Ccro  (6^4  0)1 


where  m is  the  mass  per  unit  length,  x^  is  the  center  of  qravity  of  the 

cross  section  (positive  aft  of  the  elastic  axis),  p#is  the  mass  density  of 

the  blade,  and  A is  the  cross-sectional  area  of  the  blade.  It  is  assumed  that 

the  cross  section  is  symmetrical  about  the  x axis  when  the  sectional  inte- 

P 

grals  are  carried  out.  The  position  vector  ? , the  distance  from  the  deformed 

P 

elastic  axis  to  the  point  on  the  blade  cross  section,  may  be  written  in  the 
x",  y",  z"  coordinate  system  as 

?p  « [T]  r, 

*Xp  CcrO  0 -V  Q 

~ ^p(  - u'coo  Q -V  UT'/a^Q)  — Zp  ( \r'  4 OtaO  0) 

- Xp  -V  "Zp  Cera  0 

(15) 

-►I 

The  inertial  moments  q per  unit  length  in  the  x" , y" , z"  coordinate 
system  are  given  by 

%1  - -j  J% 

A 

lixXi  I)  QC  'ir exso 0^ 4 u r 

_ ^ y*  4 r)_Au<.  0 

- 'v4  Xx  ( i}"^a-c^v0  4 UTcjcrO  Q)  - m 0„> 

-V^Xxrf  [-(r04>T)  {^0,0(0^4  0)4^000©} 

4^j  ^0^4  \Top6m 

+ ^t^M'C©m4©)J  +X0(©r4<^) 

* c T*p  ~ itp)  if  c«oCe*, + ©) 

Q.  M 'Xx  .O.  c*sO 0 C hj~  coo  0^4  vJ~ ) 

- 'Vn'X.x -rf  Cr04  v~)  coo© 

(16) 
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where  I is  the  moment  of  inertia  of  the  cross  section  per  unit  length  and 

I , I are  also  defined  as 
zp  xp 


I.  = 


= C-Xp  ■+  Zp)  d%?  dZp 


1 


■kp  = \ 


(17) 


Thus  far , inertial  loadings  are  derived  about  the  undeformed  elastic 
axis.  Next,  the  inertial  loading  about  the  feathering  axis  will  be  obtained. 

The  position  vector  r , the  distance  from  the  feathering  axis  to  the 
point  on  the  blade  cross  section,  is  given  in  the  x' ,y’ ,z'  coordinate  systems 
by 


PA 


- bp]  T Z - [Z  + [-Tp]  (Z  + [ t„]  [rs]  Y%  ) j 
= [T®]  [Ts]  ( »|  - rR  + [t]  v?  ) 


(18) 


where 


^ ^ r 


(19) 


Hence , 


20 


where  M „ and  M „ are  structural  bending  moments  and  M „ the  torsional  moment 
x z y 


V and  V are  shear  forces  and  T is  the  tension  force 
x z 


External  forces  and 
moments  of  inertial  and  aerodynamic  origin  are  described  by  p ,p  ,p  ,q  ,q 


and  q^ . Substituting  Eq.  25  into  Eq.  24  after  differentiation  of  Eq.  25  by 


r yields 

•2>3Mx» 

+ 

a r1 

ar 

a Mr- 

-V 

3 

4 

ar1 

ar 

+ ( 

'aMz 

aur 

ar 

■ 3^ 

ar 

T = 

l 

Moments  in 

Eq. 

26  will 

'd'T 


3r  J 


ar  ar y 


, BuP 
iZ  3r 


+ %X  + % 


= 0 


(26) 


formed  elastic  axis  coordinate  system  with  the  transformation  matrix  [T' ] 
obtained  from  Eq.  2 with  0 equalling  zero: 


( Mx« 

j 

i 

\j;  o 

' Ma" 1 

1=M| 

> 

i hi- 

/ 

-u 

1 -ul' 

MZ" 

V / 

0 

L 

ul'  l 

(27) 


For  the  force  components  T,  p^,  etc.,  the  same  transformation  law  holds; 


therefore,  an  identical  transformation  may  be  used  with  sufficient  accuracy. 
Applying  the  transformation  law  in  Eq.  27  to  Eq.  26  yields 


J 
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cPNU 

+ ^ | 

1 

£ 

II 

0 

3 it  ar 

< Vt-  J 

d V2 

i , 3 

3 r sr 

(y*4) 

+ Px=0 

- 1 

+ 3^  * 

+ q ^ 

4 9 — 

uX  d'T 

+ V° 

T=r  Rj<*r 


(28) 


These  are  the  usual  simple  beam  theory  equations.  It  should  be  noted 
that  the  underlined  terms  in  Eq.  28  are  elastic  torsion-bending  structural 
coupling  terms  introduced  by  Mil'  et  al.  (Ref.  41).  Similar  terms  when 
the  motions  are  rigid  pitch-flap  were  identified  by  Miller  and  Ellis  in 
Ref.  14. 

2.1.4  Moment  and  Displacement  Relations 

Based  on  the  work  in  Ref.  2,  the  moment-displacement  relation  is  given 

in  the  principal  axes  coordinate  system  (x  ,y  , and  z axes)  bv 

P P P 

- Elf  ( or"  + kr"vr°ej 

= (<;j  + £ t)  <#>'  + 4lT 

' W (29) 

= ETc  (-  vr"  crro  0 -P  uS"  a^Q)  + e*']" 


where  and  Mzp  are  bending  moments  and  Myp  is  the  torsional  moment.  Flap- 
wise  and  lagwise  bending  stiffnesses  are  EIp  and  EIc,  torsional  rigidity  is 
GJ,  tension  force  is  T,  e^  is  the  neutral  axis  location,  positive  aft  of  the 


elastic  axis  of  the  blade,  respectively,  and  k is  defined  as  follows: 

A 


A = ~~  \\  0*^)  ^Pd?f 


The  blade  cross-sectional  area  is  expressed  as  A. 

These  moments  are  transformed  to  the  undeformed  elastic  axis  coordinate 
system  (x",  y",  and  z"  coordinates)  by  the  transformation  matrix  [T] : 


= M 


,v'*P 

M zp_ 


(ETf  coo1©  -V  ETC  ur" 

-V  { ( EX^-  EIC)  c^e  A^-Q  \ -O-" 

t JA*_T  lT^r-_G_ 

+ ) c£7  + TGtw 

— 0 + VJ”  0^ 

- ( ETp  A^?-0  + ETt  c*»*-0)  or" 

~ CETy:  ~ EcTe)  C-X&  G | kT 

+ (^aT  Qtw)  uj-/  -f  e^T c*ro  G 


2.1.5  Partial  Differential  Equations  of  Motion 


By  substituting  Eqs.  14,  16  and  31  into  Eq.  26,  the  partial  dif- 
ferential equations  of  motion  can  be  obtained  in  terms  of  v,  w,  and  0 . An 

R 

additional  equation  for  the  rigid  pitch  will  be  derived  later  in  this  sub- 
section. It  should  be  noted  that  in  the  bending  equations,  only  first-order 
terms  are  retained  while  in  the  torsion  and  rigid-pitch  equations  up  to  second- 
order  terms  are  retained.  Generally,  the  third-order  terms  are  discarded  ex- 
cept when  these  third-order  terms  have  direct  effect  upon  the  torsion  and 
rigid-pitch  motions.  In  other  words,  these  are  diagonal  terms  of  the  coeffici- 
ents of  the  torsion  and  rigid-pitch  equations  if  these  equations  were  ex- 
pressed in  matrix  form.  Small  angle  assumptions  should  also  be  used  as  fol- 


lows : 

060  ©vn  = CX^©0  _ (0^-4-©^)  Q o 

©•»<.  = /3-c-n^  Q0  "t"  (©/?  "i  Oc.)  c<ro  ©o 


and 

GO"0  Q =•  CxS'O  ©Twj  — ^ Q-r*j 

Qrvl  “V  ^ ©_rW 


(33) 


where  0 , 0 , and  0 are  the  collective  pitch,  the  rigid-pitch  motion  and  the 

OR  C 

control  pitch  input,  respectively.  The  built-in  angle  of  twist  0^  and  the 
collective  pitch  angles  can  be  large.  Thus,  the  equations  may  be  written  as 
follows: 


w equation: 

^ Elp  CeS^3XG Tw  + A^-^G)  Tvo)  UT 

4 ( CEIt=.  ElJ  cooGTv>j  j 'O’"  ■+  O^T  4]’ 

“ [2  vn  Xx  Q Brw  C XT  Ccro  Qa  4 UT  Q°) 

- 'K  Xx  _Q_X  ( 4 v'")  (cuso  4]  -CTU'/ 

+ n uj  - '2'wv  .Q  ((3p  4 So  000  ^0  O' 

— vn.  ^ ^>pCv^  4y~)  /k^x  0.  4 4~ Sjj  0^  C2.0O) 

4 -lT  /vo»x©o  cjtoQo  4-  uT  /u^.2  Q0  — e0  Qk  cxro-Q  o 

— 'H  Xx  ( + ^)  cod  0T w - -\j0p5  - 2.'M.ir©(? 

“ 2 >n.ri  Xr  ^u^H0o  (-  'O  «ro  0rw  4 uj 1 a^xQt*i) 

— vw  Xx  -Q.*  | Q(5  Ctro  ( 20o  +■  Gtw_)  ~ (@o  4 Or*)  /**~i.Q0 1 

- Pw  - k (y-6  + r)  C Sj>  -V  Pf  ) 

4-  “M.  .d.  /^-Qa  (—  -V  y~ ^*D  ><3-<-v\.0o^ 

4 'M  O.  Xx  A-c-*x0o  C-crO  ^0O  4 Q-r*j) 

- | 'Wt  Xx  D.X  0TW  C ^ + O | - (eMT/4^.Grw),/ 

(34) 
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v equation  : 


C C E Ip  4 E Xc.  Cxt^QT}iJ)  v" 

c<>r°  0rw /li*<.0rM,J  Ur"  4 T C^)J 

4 £ 2 vn  Xz  il  coo  9-rw  (lr  osoQ0  kT  >0^1. 0OJ 

+ ^^t01C«rar)  (xU^Qr^J  cf>  ] , -(TuY 
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- >n.  ri  [ pyC^  + ir)  0R  ^<^G0  4 •^S'dGk  c^(20o) 

4 T_r  cuso2-0&  4 uT  /A<^_  0.  Ccso  ©,  4 <o.  /oJ^_  Qa  ^ 
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42M^nc^0o(  ir^Qn,  - ut'a^Q-tu) 

4 - CC.il  { eR  s^(260  + 0rJ  4 C/)A^-(0o-f-  0^)000.  J 

= Pv  4 ^-Q^Cn  + 'OppWG, 

4 -HU  c^o0o(  -€o  +^D/1h0o) 

4 ^a\tCOT(e64^w0t 

+ l^^ir/c^e^C^+r)}'  -(e„T A^drJ 


(35) 


26 


(t>  equation; 


[C£J  + j - (e^T /*^c0O  -(e.iTc^erw) 

+ (£i«-EJj[(Cv*)*-  (ur'')aj  A^©o  a + <*°@&)] 

+ ^^ill(r64r)  C^l'c^Qrw  4 oj'a^Otw) 

^ C'VJ"  /^Orw  -f  UT  eA°  0-rvj)  4-  "1*1  V'Sd  ©t*/  ©/? 

+ C^-ht)  PpCQr  + cI))  ^ (Qo^Qt-w) 
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- 'M  Xx  n @R  cao  Go  ( 0O  4-  0x«v  ) 
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- m oix  n ( -e„  4 s<u^.Go)  (Qr  + cro  (Q^-f-^,) 

- To  C ©^  4-  ^ ) — C^ep  ~ Txp}  TT.  C®R  ( ^-©o  + 2©rw) 
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4-  ^ Xl  T\  A«^.  (Qo  -V  Qrw)  (-e*  + ^-0o) 

+ 1 ^ Tg?  - Txp')  _QZ  a^C.^.©*  4-  2©-r»»)  — (t  -^a  0 t«j) 


(36) 


27 


where  Pw  and  Py  are  aerodynamic  forces  corresponding  to  deflections  w and  v 
(positive  the  same  directions  as  deflections)  and  is  an  aerodynamic  moment 
about  the  elastic  axis  (positive  nose  up) . 

The  quantity  T in  Eqs.  34,  35,  and  36  is  obtained  from  Eq.  28  and 
may  be  expressed  as 

T = ^ ^ <rx^(r04r-)  - 2^.0.  (xrc*£oQ04-  kS  (37) 

In  rigid-pitch  motion,  the  structural  stiffness  is  given  by  the  control 
linkage  flexibility,  represented  by  a spring  element  Kq.  The  aerodynamic 
moment  for  the  rigid  pitch  motion  is  expressed  by  Mq  (positive  nose  up) . The 
equation  may  be  expressed  as 

C ^ + ke9k  = Me  <38> 


-s: 


Substituting  Eq.  24  into  Eq.  37  yields 
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4 '*n  .O.  ^>pCv^+'r)  C0R4(gf))  C©«  + 6tw) 

->1-0.^  r(r.  + r)  &>  0R  c^G0 
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4-  **.  .O1  rSo  Xx  (zQz  + 4>)  + 
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4 V'"SD  ~ ( 20o  +©Tsi)j 
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+ A<^-Oe„)  -V  Xx«o(2Q6+  QTtg)  J 

4 ccrt(2Go)  - -^-  'U'1/A^_(20o)  4.  -4  ovvQ^1A^(^0')Jrl/' 

^ ^ £"rt  C^P  ~ 1*p)  Q.  A4'*'-  (200  + ^Qxw) 

-\-  /vnllapp  (v2  4r)Q:xaso(eo-VeTw) 

4 ^p  v-C  rc  + 'r)^T>  - 'H.Cfe.  v^S©  0,50  ^ 

-V  "wv. Cl?  Go'Xx  /x-^v-CGo-f  0 Tw) 
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(39) 


The  final  equations  (Eqs.  34,  35,  36,  37,  and  39)  are  the  basic  equations 
of  the  nonuniform  rotating  beam  and  are  nonlinear  partial  differential  equa- 
tions with  nonconservative  (aerodynamic  and  mechanical  damping)  forces. 
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2.2  Natural  Frequencies  and  Associated  Mode  Shapes 


2.2.1  Introduction 

In  this  section  the  finite-element-method  formulation  to  obtain  the 
natural  frequencies  and  associated  mode  shapes  is  described  briefly.  The  formu 
lation  itself  is  rather  straightforward:  the  coupled  motions  can  be  expressed 
in  terms  of  the  natural  frequencies  and  associated  mode  shapes.  It  is  expected 
that  coupled  natural  frequencies  and  mode  shapes  will  give  much  faster  conver- 
gence than  uncoupled  standard  mode  shapes  for  a nonrotating  uniform  cantilever 
beam  will  do.  The  faster  convergence  will  be  pronounced  when  the  motion  be- 
comes strongly  coupled  between  flapping  and  torsion  due  to  blade  chordwise 
mass  balance  or  between  rigid-pitch  motion  and  elastic  torsion  due  to  control 
linkage  flexibility. 

Although  it  has  been  recognized  that  the  finite-element  method  is  a gen- 
eralized Ritz  method  (Ref.  42) , the  finite-element  method  is  much  more  versa- 
tile because  in  the  Ritz  method  the  assumed  displacement  modes  are  extended  to 
the  entire  structure;  thus,  for  the  structure  with  several  different  boundary 
geometries,  it  will  be  difficult  to  choose  appropriate  admissible  displacement 
modes.  In  addition  to  this,  when  an  accurate  solution  is  required,  a large 
number  of  assumed  modes  must  be  used  and  considerable  algebraic  manipulation 
must  be  performed  in  the  Ritz  method.  In  the  case  of  the  finite-element  method 
the  same  assumed  displacement  function, which  is  usually  simple, can  be  used  for 
each  individual  element,  and  boundary  conditions  of  several  different  types  can 
be  applied  simultaneously  without  difficulty. 

To  obtain  the  natural  frequencies,  the  second-order  differential  equa- 
tions without  aerodynamic  forces  are  solved  which  consist  of  the  symmetric 
mass  and  stiffness  matrices  and  the  antisymmetric  gyroscopic  matrix  due  to 
coriolis  forces.  This  problem  is  an  eigenvalue  problem;  however,  mode  shapes 
must  be  expressed  by  complex  numbers  expressing  the  phase  lag  between  motions. 
This  creates  further  complexities  in  the  complex  number  algebraic  manipulation 
for  a Galerkin  method  application.  Also,  since  the  coriolis  force  effect  is 
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not  strong  enough  to  make  a difference  in  convergence  between  cases,  it  is 
included  in  the  natural  frequency  and  mode-shape  calculation,  excluded  in 
the  natural  frequency  calculation,  and  then  included  later  in  the  Galerkin 
method  application  step.  Therefore,  in  the  present  analysis  the  natural  fre- 
quencies and  mode  shapes  are  obtained  from  the  eigenvalue  problem  with  only 
the  mass  and  stiffness  matrices,  which  will  produce  real  mode  shapes  and 
natural  frequencies. 

2,2.2  Variational  Functional 

For  the  first  step,  the  variational  functional  it  will  be  derived  from 
Eqs . 34,  35,  36,  37,  and  39.  Nonlinear  terms  and  steady-state  loading  terms 
are  excluded  as  well  as  velocity-dependent  terms  resulting  from  Coriolis 
forces.  The  potential  energy  and  kinetic  energy  expressions  related  to  the 
linear  terms  are  retained: 
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-kOz  (pP  C Y4  4 'r~)  c-<f=©0  4 V~  £Tc>  c*ro  (^2.0o) 

-<-e«^~-©0  - Xx4m^U©.4  0tJ  j *ir^ 

+ {^'Xx  AX r~So  ^(20„4B^) 

~ ^ ClXQ0  ^ C©o4  6-r  J 

4 (t2p  -T  xjp')  Q?  wra  ( 2&0  -v  2.0Ti(j) 

2^  ^<©  0^  (40) 
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2.2.3  Finite-Element  Formulation 

In  the  finite-element  formulation,  the  variational  functional  may  be 
expressed  as  the  summation  of  the  functionals  of  the  subdivided  finite  number 
of  discrete  elements. 

N 

TT  - J"  7T^  (4i) 

/M=l 

where  n means  the  nth  element  and  N is  the  total  number  of  elements.  The 

displacements  |uj  are  represented  over  the  nth  element  by  interpolation 

functions  [f]  and  by  generalized  nodal  displacements  {q  } which  are  dis- 

n 

placements  and  sometimes  derivatives  of  the  displacements  of  the  nth  ele- 
ment. Thus,  in  the  matrix  form 

M = [i  ] I <42' 

where 


L u = u.  UT  xT~  4>  0R_, 


;+l  tj;  tj-khi  nj'X 


%m  -i  t_ 


(43) 
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In  the  present  analysis  the  beam  element  is  used,  and  w and  w'  are  the 

n n 

vertical  displacement  and  its  derivative  at  the  node  between  (n  - l)th  ele- 
ment and  nth  element.  The  quantities,  w , and  w'  , are  those  between  the 

n+1  n+1 

nth  and  (n+l)th  elements  (Fig.  6).  The  matrix  f is  written  as 


u(" 

M0I 

h!? 

u<“ 

n01 

W2 

0 

O 

0 

o 

0 

o 

o 

o 

O 

0 

o 

0 

0 

HJ 

Ho’* 

m2 

o 

o 

o 

o 

O 

0 

0 

0 

0 

O 

o 

0 

0 

u°' 

“oi 

h*;? 

Ho  ‘l 

H?i 

o 

0 

0 

0 

O 

0 

o 

o 

o 

O 

o 

o 

o 

1 

where  the  Hermitian  interpolation  functions  are  defined  as 


(44) 


H2  - >(i)3  - 3(1)%  I 
< . -*(*)*  4* 
h‘,?  = 


The  element  size  is  expressed  by  Z,  and  x is  the  coordinate  within  the 
element.  Two  nodes  of  the  element  are  located  at  x = 0 and  x = Z. 

The  derivatives  can  be  obtained  from  Eq.  42.  For  example, 

aur- 


up  + ^ tin  uy/  ■ d^'i 

dX  " dZ  ”M  dx  ^ + 


dx  ^ 


(46) 
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Thus,  various  derivatives  can  be  obtained  by  differentiating  interpolation 
functions.  These  displacements  and  derivatives  (Eqs.  42  and  44)  may  be 
substituted  into  the  functional  in  Eq.  40  and  the  integration  should  be 
carried  out  elementwise. 

where  [m  1 and  [k  1 are  element  mass  and  stiffness  matrices, 
n n 

The  element  mass  matrix  is  a consistent  mass  matrix  and  may  be  written 


d„0']  to] 

s [Wj 

13* 

symmetrical 


d,3  ['*<] 

da  Y.™'} 

dU-,  t 

da  W'J 

oLh  £ 

o-j = J_ 

J -fao 


156  12(1  54  - 13d 


■ rvi-I  e 

is&  -2i(  ’ L J 12  6 


symme trical 


The  coefficients  are 

^ ~ , C^I3=  - ItC-OgQtw  , C*rO  Q-th, 

^ -‘M'U  (50) 

C^33  "^°  ; <^3-v  = Xo  - 'Hi  Z.x.  /z^nQ-rw 

d^4  as.  To  + — Q 'Hi.  X-z  <^D  /0-Ly*i.Q TVJ 

r : radial  distance  from  the  pitch  bearing  to  the  midpoint  of  the  element. 
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The  element  stiffness  matrix  is  divided  into  three  portions  as 

[«j  = [<:]  + [*:]  *[*i] 

The  stiffness  related  to  the  curvature  is 


IX]  = 


*"IX]  sia[4;]  Co]  Co] 

s^[4s]  [0J  [Q-j 

symmetrical 


W]  = 4: 


\i  £$  -IQ  Qi 
-C£  2^ 


symme trical 


IQ 


Sh  = EXp  c<W10Tw  ■+■  f lao^'Qrw 

sii  = CEtp  - EIc]c«r3  0TK/<»j^i.0Ttv 
Sii=  EXp  ^w.l©rvj  + ETd  OSO^Gtiv 


The  stiffness  due  to  centrifugal  forces  and  torsional  rigidity  is 

^r4]  to]  Lo]  Co] 

IX]  - tJ2Ut]  Cp]  Co] 

Co] 

symme trical 

L O 
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3 Q 3$  -a6  3^ 

4 1 -r 

sQ  — 3$ 
symmetrical 


t„  = ~T 

rR 

= l >M.C.<;+r  )£1 

^ U 

"^22.  = T 
tjB  = ^vJ  4 


The  stiffness  due  to  centrifugal  forces  and  mass  unbalance  is 


[*"]  = 


a„  CVJ  aa  [">*.' ] 

O] 

(3-aa  ^ '*'1  J 

ai2>  [Wj 

a«L*<] 

(56) 

t0*1'] 

a*S  [ />4"J 

symmetrical 

0U<v£ 

Qo  Q. 


a"  = - 'K-O.'1  a^2Q0  , au=  ^.Qaat to06 

a,3~  ^Xia^CG^  + G^^e, 
l<V  ~ ^ O | pp  Cvb  + |r^j  -+  rK  £Tj>  ,4^(20.)  - Qo 

+ 4 QTly)  j 

VtO.  c^=  0O  , <lii  = ^ 'Xx -O!1  CGo+  ©Tw)  C^dQ0 
{ Pp  C V.  + r-.l  ^e„  4 ^ ^g.)  4 e.^G. 

~ “Xx  C 2©»  "V  ^ 

~ Cl*p  'Ixp)  <^°  ( 2-Go  + 2.Qtw> 

4 ^C^  + Grw)  (~e.  4 KT<^/a4t_^0) 


v*'  Qtw)  4 MXi  O.1  C*ri+v'0 

Qy*  = ( I^P  - X^_p3  Xi7- cur=3(2©e,+  2©rw)  + ”M  Ii  /^-(20.1-Sw) 

— >*\  lx  .£"1*  e 0 c jsv(jQ0  +©rw)  — >n  Xx  SnlC^  ’r*.)  /0^^-C®0^  Q*w) 

^44  =■  ( Izp  — Xxp) -O  C<ro  (2©o  + 20rw) 

- im.  .Q*  | (3p  v'^.)  c^°©«  4"  2 €?0  /''»  c)c> 

1"  Q.  &*  e^=  ( 2©„)  4 2eoX1Mt>(Q64  0T«) 

— *4  X'h  Sp  ^ x ( 20*4  0ru)  | 

J (57) 


Substituting  Eq.  47  into  Eq.  41  and  assembling  the  element  mass  and 
stiffness  matrices  appropriately  for  compatibility  of  the  interelement  nodal 
displacements,  the  global  mass  and  stiffness  matrices  are  obtained  and  the 
expression  for  T r can  be  written  as 

TT  -7  {?fMUi  - 

where  {q}  is  the  global  nodal  displacements,  and  [M]  and  [k]  are,  respectively, 
global  mass  and  stiffness  matrices.  These  may  be  written  as 

Ite}  Kj {«.}]-  WUI 

[mJ  = [j]  [[•«.,]  [h,] [-W.J]  [j]  (59> 

[k]  = [o\f  [[4.H4J  [*„]]  [j] 

Matrix  J is  a transformation  matrix  to  relate  the  element  nodal  displacements 
{qn)  the  element  to  a column  of  global  displacements  {q}. 
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The  application  of  Hamilton's  principle  to  Eq,  58  will  yield  the  equa- 
tions of  motion, 

•+  Ik]  (si  = {°J  «6(>» 

2.2.4  Boundary  Conditions 

Various  boundary  conditions  can  be  chosen  depending  upon  the  rotor  con- 
figuration and  can  be  applied  easily  in  the  finite-element  method  to  specify 
conditions  of  the  displacements  of  the  first  node. 

Hingeless  Rotor: 


-UT,  = Ur/  _ Tj;  _ qj;'  _ o (61) 

Gimballed  Rotor  (Teetering  Rotor) : 

Uf,  = \T,  = ^ = 0 (62) 

Articulated  Rotor: 

^ = OT,  = _ o (63) 

If  control  linkage  flexibility  exists,  the  spring  and  constant  Kq  should 
be  included  in  the  element  stiffness  matrix  of  the  first  element.  If  the  con- 
trol linkage  is  rigid  enough,  the  rigid  pitch  degree  of  freedom  should  be  re- 
moved before  the  eigenvalue  calculation. 

2.2.5  Subspace  Iteration  Method 

The  eigenvalue  problem  expressed  in  Eq.  60  can  be  solved  by  any  eigen- 
value technique.  In  the  present  analysis,  the  subspace  iteration  method  (Refs. 
43  and  44)  is  selected  to  achieve  a fast  rate  of  convergence  of  a large  degree- 
of-freedom  system  (if  ten  elements  are  used  to  describe  the  rotor  blade,  the 
total  degrees-of-freedom  will  be  67) . 

The  subspace  iteration  method  will  be  described  briefly.  Assuming  s 
modes  are  needed,  the  total  degrees-of-freedom  as  t ({q}  in  Eq.  60  is  a 
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f 


t x 1 matrix),  and  denoting  [u]  as  the  t x s matrix  containing  the  s vectors, 
the  following  procedure  will  be  applied: 


a.  Assume  an  arbitrary  [u  ] as  a starting  matrix 

o 

b.  Calculate  the  reduced  mass  and  stiffness  matrices  as 


[mJ  = 

uf 

[m]  [Uo] 

5XS 

xx  t 

t*t  *”s 

and 

[kr]  - 

[-r 

Lk]  LwJ 

SxS 

sxt 

txt  T*S 

Solve  the  small 

eigenvalue  problem 

Lk*]  La]  - Ld][. La]  - (o } 

where  [A]  is  a sxs  eigenvector  matrix  and  [D]  is  a 
diagonal  matrix  which  includes  eigenvalues.  This 
small  eigenvalue  problem  can  be  achieved  by  using 
the  Jacobi  Method. 


d. 

Calculate 

^2 

13 

1 1 

as 

L*o] 

= [0 

[A] 

t*  s 

SXS 

e. 

Calculate 

as 

[u,] 

= [k]-' 

[ M]  [ 

f . 

Replace 

Lu°] 

with 

[u] 

and  go  to 

step  b. 

The  criteria  of  terminating  the  iteration  is  set  to  achieve  the  specified 
accuracy  as 

0);.,  - O); 

where  (AK,  (AK+^  are  ith  and  (i+l)th  iteration  eigenvalue  and  e is  a speci- 
fied accuracy. 


< £ 


(64) 
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2.3  Aerodynamic  Loading 

2.3.1  Introduction 

The  aerodynamic  forces  and  moments  of  the  rotor  blade  are  derived  based 
on  the  two-dimensional  quasi-steady  aerodynamics  in  Ref.  45.  The  major  assump- 
tions in  the  present  aerodynamic  analysis  are: 

(1)  the  inflow  is  uniform  over  the  rotor  disc  plane,  and 

(2)  the  vertical  gust  is  assumed  to  be  one-dimensional  in  cruising 
direction, 

(3)  reverse  flow  and  blade  stall  are  neglected:  thus,  the  analysis 
is  good  to  an  advance  ratio  of  about  0.4. 

Also,  the  rotor  unsteady  wake  effects  are  neglected.  The  virtual  mass 
aerodynamic  forces  and  moments  are  neglected;  the  terms  of  order  c/R  in  the 
aerodynamic  lift  expressions  are  neglected;  the  firstorder  velocity  terms  in 

3 

the  aerodynamic  loading  are  retained;  terms  of  order  (c/R)  in  the  aerodynamic 

feathering  moments  are  also  neglected;  the  collective  pitch  0 and  the  built-in 

o 

angle  of  twist  0 are,  hereafter,  restricted  to  small  angles  which  enable  the 

xW 

use  of  the  small  angle  assumptions  of  helicopter  aerodynamics. 

2.3.2  Relative  Velocity  of  the  Blade 

The  position  vector  r to  the  point  on  the  blade  elastic  axis  in  the  x,y, 

B 

z coordinate  systems  after  applying  the  small  angle  assumption  as  0 <<  1 and 

m 

0 <<  1 to  Eq.  11  is  expressed  as 


-e„  + V 
rv 


UJ-  4-  T'^D  “VC  v-o'+'O 


(65) 


41 


The  velocity  due  to  deformations  may  be  written  in  the  x,y,z  coordinate 
systems  as 


-v  <fi  * 

■\T  - .O.  "V- f 0 ( + y~S\3 ) 

« 

nc-e.  + ir)  + Q ( -e0+  or) 
- e(-e„  + -Ur) 


(66) 


where 


0 


0 

•Q.+ 


(67) 


and  0 includes  the  angula.  velocity  due  to  the  rigid  pitch  motion  and  the 
cyclic  pitch  control,  while  the  blade  elastic  torsion  has  no  effect  on  the 
velocity  at  the  blade  elastic  axis  because  it  is  assumed  that  the  elastic 
torsion  deformation  occurs  after  the  vertical  and  chordwise  bending  deforma- 
tion shown  in  Fig.  3. 


The  velocity  seen  on  the  projection  of  the  blade  elastic  axis  to  the  x-y 
plane  is 


= orQ  - Cl  x ^ 


\r  - XlS  -Vvi  d-  e C vr  4-  ) 

C -Q-  + 0)  \r 

- ur  - e ir 


(68) 


The  trim  velocity  may  be  expressed  in  the  x,y,z  coordinate  system  (the 
hub  plane)  as 


42 


V = 


|LA.  XI  R 

p'  XI  R CJTO  (p 

- a xiR 


where  p,  X,  and  ag  are  the  rotor  advance  ratio,  inflow  ratio  and  the  rotor 
shaft  tilt  angle  (positive  rotor  disc  tilt  forward) , respectively,  and  de- 
fined as  follows: 

a V O CTO  0 (s 
r ~ xirR 


A = 


~\f  Xs*  “I"  '"'A 
XI  R 


The  cruising  velocity  is  V,  and  the  induced  flow  velocity  is  v . The  resultant 

, . - 1 

velocity  U may  be  written  as 


U = V - 


+ fieri -(-r)  - vr~ 8(io-  + 7f+r5) 
^Q.R  ocrO^>  - XI  ( 

-A  XIR  - UJ-  -0  (xr  + X?) 


This  resultant  velocity  U is  resolved  to  obtain  the  relative  velocity 
components  of  the  deformed  elastic  axis  coordinate  systems  (x''',  y' ' ' , z' ' ' 
coordinate  systems  shown  in  Fig.  3)  as 


u, 


£4 


U* 


"Up 


= wT[TdT[iPru 


jO.Cy'0+v^’')  4-  f^QR  ~ 1-r 

— p'  O-  R XT1  curO  Cp  - 0Cul4-?p+^) 


-XQ.R(ul  + p>p+^) 
+ f^Q.R  ar,/^t-c{J  -4-  ciC^.-vR)  u-' 


- - Rr  -p.£XR(uj-'  + (3^+^)^°}' 

+ 0C4>  + Xf) 


(72) 


where  U , U , and  U are  blade  tangential,  radial,  and  perpendicular  veloci- 


T 'R'  p 

ties  with  respect  to  the  deformed  elastic  axis. 


The  bending  motions  are  expressed  in  terms  of  steady  equilibrium 
quantities  and  small  unsteady  perturbation  quantities 


or  = 4-  \r 

ur  vTc  -+  ur 

The  same  symbols  are  used  for  the  total  bending  motions  and  perturba- 


tion motions;  however,  there  will  be  no  confusion. 
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Equilibrium  velocities  related  to  the  trim  condition  are  obtained  as 
follows : 

r U-, 


U 


iJ-R. 

V-Up  j 


- qj;  -HClRvz'c* ^ 

ixaRc^oc/.  -ci(^+*p)- 

+ 'u;/xuUt_c/)  + JQ.C^ +r-j 'U-o' 

-AQ.R  - uTs  + 


(73) 


where  v and  w are  bending  deflections  at  the  equilibrium  condition  shown 
o o 

previously,  and  time  derivatives  v and  w are  not  necessarily  zero  because 

o o 

in  the  present  analysis  no  trim  operation  is  considered  which  involves  time- 
dependent  equilibrium  position-related  harmonic  terms;  i.e.,  siniji  and  cos \(j. 

The  perturbations  of  the  velocity  components  including  the  vertical 
gust  are: 


six  = 


% 

- ir  - ^ HR  'it1  cjoo<4>  — 9(  ur+ ^&t>) 

■+■  /O^tCXa  0 

- Qir  - A nR  us'  4-  /uQR 

_ _ + jSp.  C_r£  +irj  TJ- ' -\r_  UTS  C^°_CP  _ 

~ US  — [A,  -O.  ^ uj'  euro  (J)  +■  @ tT  -f  ^ 


(74) 


where  w is  the  vertical  gust  at  a given  location,  positive  upward. 
G 


2.3.3  Section  Aerodynamic  Forces 

The  section  lift  and  drag  are  (Fig.  8) 

L = 4:  J3  o.o( 

D = 4 c*o 


(75) 


where 


-LT  = UPl 

o = e - 4>x 


(76) 


Cf>t  = -W1  -S.  „ 

UT  U-r 


and  a is  the  lift  curve  slope  and  the  profile  drag  coefficient.  The 


force  components  P , parallel  to  U , and  P , parallel  to  U are 
z p x T 


Pz  = L - D + 4 pa(§f  UT  (3+4TCA)(?w) 

=*  -Lp<kc  UUTo(  -^-pca0  UUF 


+ ut  o + 4**k£°<) 

+ D c^<px  +-4p°.(-%f  UpCa-v^^OC &*) 

— U p a c TJ  Up  & -+  4 P c C,k' 

4-lpcx(%f  XJp  (3+4^  C&x) 


(77) 


where  xft  is  the  nondimensionalized  distance  of  the  blade  chord  between  the 


elastic  axis  and  the  aerodynamic  center  (positive  aft  of  the  elastic  axis) . 

It  should  be  noted  that  the  aerodynamic  force  in  the  blade  radial  direction  is 
neglected  and  the  lift  due  to  the  blade  pitching  angular  velocity  is  included 
to  retain  the  damping  terms  of  the  torsion  and  rigid  pitch  motion  shown  in  a 
later  section. 
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The  transformation  [Tq]  in  Eq.  6 may  be  used  to  resolve  the  blade  forces 
P2  anc*  px  into  Pw  and  parallel  to  the  z"  and  x"  axes  (the  directions  of  the 
bending)  of  the  undeformed  blade  coordinate  system. 

= P2  + Px  Qm 

(78) 

P„  - - B 6m.  + Px 

The  blade  section  velocities  and  angle  of  attack  are  expanded  into 
equilibrium  terms  and  perturbation  terms  as  follows: 

UT  — UL-y  4-  £uT 


Up  = 


up  4 <^Up 


XJ  = -4  g'u-r 


= C.  0o  + 0rw)  4^0+  Sd>  - — SuT 

“ u*-ru* 


where  uT»  up,  6uT,  and  6u^  are  defined  in  Eq.  73  and  74,  and 

U = Jup  + ul 

Se=  Gk  -h  0c  ; = c (80) 

Then,  the  aerodynamic  forces  are  resolved  into  equilibrium  terms  related  to 
trim  operation,  and  perturbation  terms  due  to  perturbed  motions  of  the  system. 
Substituting  Eq.  79  into  Eq.  78  with  Eq.  77  yields 


P»  - <*-)-■*•  (f§0.^  +(l|)„K+(#).(S‘e+Sct=)-(§!)>Csre+S^) 

- <R).  + (ftj.  1H  +(4|)„ » + 

+ (f|L  («  + *P) 


where 

(Pw)o  = -1  pcx C (£1R T [ A20s  4 xpx0sx^4 

*t  p*  0'S  f)  — A p A*^-^  ~ C \ + ) X (A  4 ^)j 

(~f~)  = 4;  PaC  C-Q-R)  F^TcQs  4 Ap0s  '<^4 

— M-  (uJ'.l+  ^j>  -V  St>)  c«o</>  — ( \ 4 (>\4  Sro 

*-2P°-c(£iR}[0  + a.)  Cx+H^A - "Oj 

~ A <^TO  Cp  )^J 

(l^)&  =4  PO-cC^lR)1  [ xz  *p'X^</’ 

~ ‘i  P*  CisoQtp)  _ 2 XU;  4 2./LX  X CJS^cjj 

-IfJlU  A^</>  J 

= 4.  P aC  (r)  (4  + '*■*)  (&-£■)  [^X  + p tjj 

- ~ o^4J 

^ v^°  ~ P ^ C C-^-R)  f"  X 0s  C A ■*"  ^ ) 4 £4sl  x1"  — Ax 

L £X 

- 2XXSI  ~ 0b  | X203  4 ip  X0S  ^cfj 

+ p-1 0-s  <A*^a'4  - X /A  - C A 4 uj;  ) j j 

^rl  2 P^C^R)  [ 03  [a  4 kj;  -t-p.  ( UJ^'  4 + Sp)c<^>cj)  , 

4 X^S  X-©6{2'*©s42-PGs'a~4  “(A4*j)| 
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(■|^)  = ~ pac  (Q.R)  xQs  + p©&  >^4  - 2 1 a + id: 

-+  p. c uj*  + pf-f  g^j  oso4) j + a,  C a -f - p <^4  _ u ~p 

("§©i  = 4 P^*2  CUR)1  [ a * -+-  a:  uj. 

+ f*  "*  C_kT0/  + (^p  + S'.)  curoC j)  + 4 p2"  (" ul^  4^>p+  Si>)  /^-('24j) 

~ | Qs  + 2p  XQs  ,4-<— <-t|>  — X ( A 4-  &;) 

+ p1  ©s  - A p ^K_c^  J - Go  { A1  4-  Y P1 

-t-  ^pA.  A^.c\* i^c«ro^2c|j)  - 2/X/U; 

4-  2-p  A trj  osoc^  - ip  u;  /AUvCpJ  J 

)0  = *2  P a-<~  CUQ.)  [a'J;  4-  A id;  -t-  Ap  £u*i_<\>  + p u_T  a^-^P 

+ px(ujy  +(3,p-v  s,>)  oso  c|>  + 4 P ( 4-^p^Sfc)  A***-  (2-4) 

- ©„{  *.*■  + 4 pl  ( 2 p 2 ^^.4  - 4 p1  c«so(2.4)  - 2A  IT. 

4-  2 p x u/ ego <4  _ 2p  15;  4j] 

(Hi  =4p°-c©  (.J  + pjCaR.)  [*  * A- 

+ p ( 4-  4 5d  ) acrO  cfj  J 


It  should  be  noted  that  the  double  underlined  terms  are  due  to  deforma- 
tions at  the  equilibrium  condition.  The  nondimensional  radial  coordinate  x 
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(83) 


and  angle  ©s  are  given  by 


x = O^-^Vr 

03  = 0o  ■+■  0-rvj 


2.3.4  Aerodynamic  Moments 

The  section  aerodynamic  moment  (positive  nose  up)  may  be  written  as 


MA  = > 


P a U7e  TU  CX  4 4:  PXJ2  C1(ClHo4  c^ot) 

-4  P0,u  c3  cMcx 


where  Cm  and  Cm  are  section  aerodynamic  moment  coefficient  and  section 
o a 

dynamic  moment  slope,  respectively.  C is  a function  of  x : 

M A 

r'  1 *3  _ x 

^ -v  "Xa  + "Xa  (f 

Expanding  as  before  in  equilibrium  and  perturbation  terms  yields 

MA . (mj.  + css^sut  + (IS-;)  Sup 

+ (^](Se+S<t>) +(!§»]  (sSV**.)  « 


where 


( ma)o  = 4 pa c1  (AR)1  [-L  Am ^4 

~zpl  <^>(2.cp)  _ 2.^]  + •e>m  x.  (a+  kJ^)J 

= T pGLC  [ AnCx  + _ xpj 

_•&*(*  + Crej] 

taup)e  = 4Pa^C^R)[AM  Ca+wc) 

4 C X 4 p 


and 


2>Ha\ 
2>©  /< 


pGL  C1  ( iTiR)  ^ | 'X1  + "5?  W + 2_p.  X <<*■*—  cjj 

— T ^ ^ C^)  ~ 2^.  jj 


^3fl  /c 


= - 4:PaC-3  [cH(  - U;)J 


= 2.  (%■+%*  0s  - ^Qs) 


(87) 


?>M  — TU  ~ 

cx 


(88) 


Aerodynamic  moments  Mg  for  the  rigid-pitch  motion  consist  of  two 

R 

portions:  one  due  to  the  section  aerodynamic  moment  derived  previously  and 
the  other  due  to  lift  and  drag  which  make  moments  about  the  blade  feathering 
axis : 


1 ! 


2.3.5  Gust  Velocity  Gradient  Due  to  the  One-Dimensional  Gust 

As  shown  in  Fig.  9,  the  one-dimensional  sinusoidal  gust  in  the  cruising 
direction  has  a velocity  gradient.  The  gust  magnitude  at  the  nondimensional 
spanwise  location  x of  the  blade  of  azimuth  angle  ij;  = fit  is  expressed  as: 
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(93) 


— U3^-r  C ) 


where  w is  a gust  amplitude  and  a)  the  gust  frequency,  <J>  phase  lag  of  the 
GT  G 

gust  between  that  location  and  the  center  of  the  rotor.  The  phase  lag  <j>G 

is  described  by  the  gust  wave  length  X and  the  location  of  the  given  point 

G 


4>«  = 27T  [- 


Tc.  CUTO  C 


= -t 


Substituting  Eq.  94  into  Eq.  93  and  expanding  the  trigonometric  function 


yields: 


( CbO  oft  ca o at 


+ Higher  Order  Terms  ( 

In  order  to  obtain  the  analytical  excitation  function,  the  approximate 
Eq.  95  is  used  instead  of  Eq.  93.  The  first  term  shows  the  uniform  part 
of  the  gust  and  the  second  the  gust  velocity  gradient  due  to  the  one- 
dimensional sinusoidal  vertical  gust. 
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2.4  Rotor  Trim  and  Steady-State  Deflections 


2.4.1  Introduction 

In  this  subsection,  an  approximate  method  to  obtain  the  rotor  trim  con- 
dition and  steady-state  blade  deflections  will  be  described.  These  are  re- 
quired to  evaluate  the  aerodynamic  coefficients  and  structural  couplings  of 
the  perturbation  equations.  There  are  two  rotor  trim  conditions:  one  is 
powered  trim.  This  simulates  the  actual  helicopter  trim  and  requires  force 
equilibrium  in  vertical  and  horizontal  directions  as  well  as  moment  equilibri- 
um in  pitch  and  roll.  The  other  is  moment  trim.  This  simulates  wind  tunnel 
testing  of  the  rotor  where  moment  equilibrium  in  pitch  and  roll  is  required 
and  the  thrust  and  drag  are  reacted  by  the  wind  tunnel  model  support. 

In  the  present  analysis,  both  moment  trim  and  out-of-trim  operation 
are  considered  with  respect  to  the  wind  tunnel  experiments.  In  out-of-trim 
operation,  the  blade  motion  is  periodic.  In  this  case,  the  blade  torsion  and 
rigid  pitch  motions  are  neglected  when  the  steady-state  trim  condition  is 
obtained. 


2.4.2  Equations  of  Motion  for  the  Steady  Equilibrium  Condition 

Discarding  terms  related  to  the  elastic  torsion  <t>  and  rigid  pitch 

motion  0 , and  second-order  terms,  two  sets  of  equations  are  obtained  from 
R 

Eqs.  3 and  4 as  follows: 
w^  Equation 

= (96) 
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Vq  Equation 

i'Vh.  4 *2.  'M.  O Cpp+  Sj>)  UT&  4 -4  'vm.-Q.  'Xx  'Ol 
- "M.  --w^O1  u£ 

£ ({=1?  A^Qtn-,  + £"I^ccoi0rw)  TJ;  + {(Elfr-E^aeoQru  /•^-0Tk^  USl  J 
= Pv  4-  'M.  -R2  C vto4  |3f0G  4 -Vvt  03  ( — e„  4 <*§^©o) 

4 2.  kv>.  n.1  — 'V*t  ClzeM 

(97) 


Galerkin’s  Method  is  applied  to  Eqs.  96  and  97  to  obtain  the  modal 
equations  of  motion  for  the  steady-equilibrium  condition,  using  the  two  coupled 
rotating  bending  mode  shapes  obtained  in  Subsection  2.2.  The  dimensional 
bending  is  expressed  in  terms  of  a series  of  generalized  coordinates  and 
dimensionless  mode  shapes  which  have  unit  displacement  at  the  blade  tip: 


*01 


= R 


3.  f ) _ 

> j ' ^Ct) 

V*  (x) 


-6.-1 


(98) 


After  substituting  Eq.  98  into  Eqs.  96  and  97,  the  application  of 
Galerkin's  method  with  coupled  mode  shapes  is  as  follows: 


Left-hand  side  of  Eq.  96  in  terms  of 
generalized  coordinates  and  mode  shapes 


Left-hand  side  of  Eq.  97  in  terms  of 
generalized  coordinates  and  mode  shapes 


dr 


Right-hand  side  of  Eq.  96  in  terms  of 
generalized  coordinates  and  mode  shapes 


Right-hand  side  of  Eq.  97  in  terms  of 
generalized  coordinates  and  mode  shapes 


dr 


(99) 
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It  should  be  noted  that  aerodynamic  forces  P and  P in  Eqs.  96  and 

W V 

97  are  evaluated  from  Eqs.  74  and  82  discarding  terms  related  to  torsion 
<J)  and  rigid-pitch  motion  0 and  double  underlined  terms  which  express  the 
steady-state  equilibrium  deflection  effect  on  aerodynamic  forces;  they  will 
be  included  in  the  perturbation  equations  later. 

After  some  algebraic  manipulation,  the  equations  of  motion  for  the 
steady-state  equilibrium  will  be  obtained  in  terms  of  generalized  coordinates 
as  follows: 

+ (2S«,  + (2su, + Mu.)"fi 

'V  = Mc.j,  4 Mti  + e 


+ (2.St4i  4 i 4 (2%j,+  *1*') 

+ % Qsrf 


(100) 


where  qtand  q^  are  generalized  coordinates  corresponding  to  the  lowest  fre- 
quency flapping  and  lagging  motion,  and  (°°)  and  (°)  are  nondimensional  time 
2 2 

derivatives  d /d(fit)  and  d/d  (fit).  Iq^  and  Iq2  are  generalized  masses  and 
Sq^q^ , Sq2,ql  and  ^2^2  are  related  to  Coriolis  forces.  Rotating- 

blade  natural  frequencies  are  expressed  by  to  /fi  and  to  /fi.  Mq  q , Mq  q , 

• • J-  ^ 1 1 1 4 

^2^1  and  Mq2^2  are  aerodynamic  dampings,  and  Mc  F q and  Mc  f q are 

loading  due  to  centrifugal  forces.  Mq^  and  Mq^  are  aerodynamic  loading,  and 
M „ and  M Q are  coefficients  of  blade  pitch  control.  These  coefficients 

qr  q2® 

are  defined  in  Appendix  A. 


r ■ 


Generalized  coordinates 

q-L  and  q2  may 

be  written  as 

(30 

(3,t  -v 

(3rs 

A— 

(101) 

%(*>  - 

So  + 

Sit  ^4  -+ 

where  it  is  assumed  that  the  flap  motion  is  predominant  in  q^  (t)  motion  and 

the  lag  motion  is  predominant  in  q„  (t)  motion.  8 , 8,  , 3 , C F,  , and  F 

2 o lc  Is  o lc  Is 

are  constant  and  express  the  amplitudes  of  the  collective  and  cyclic  motions 
in  the  steady  equilibrium  condition.  Trim  pitch  is 

0 = 0o  + ©.c  + Qi*  (102) 

where  9,0,,  and  9,  are  constant, 
o lc  Is 

Substituting  Eq.  101  and  the  time  derivatives  of  Eq.  101  and  also 
Eq.  102  into  Eq.  100  and  applying  the  harmonic  balance  method  yields  the 
following  linear  simultaneous  equations  for  the  out-of-trim  operation. 


FOR  NO  TRIM  OPERATION 


(103) 

where  coefficients  used  above  are  defined  in  Appendix  A. 

For  the  moment  trim  operation,  to  achieve  the  state  of  zero  pitching  and 
rolling  moments,  the  following  condition  is  necessary: 
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(104) 


^>ic  =•  “ O 

which  requires  trim  cyclic  pitch.  Therefore,  the  equations  for  the  moment 


trim  operation  are 


( \ 
ft 


G, 


KJ 


6lS 


5J 


o 


M 


M 0 T-t, 


o 


M, 


(105) 


where  the  coefficients  used  above  are  defined  in  Appendix  A. 
2.4.3  Thrust  Coefficient 


The  thrust  coefficient  CT  is  defined  as 


C-r  = 


PCTTR’-XxiR)2 

— 5 \ \ 4)  dr dd 

^KaQ)1  )0  )o  T 


2'nja(TTR*)(SlG)7 


(106) 
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where  T is  thrust  and  N is  the  number  of  blades.  Assuming  P and 

substituting  Eqs.  74  and  82  into  Eq.  96  yields 


cr-a. 

2. Xe  ^ - Oh- -%*)  A -+  p-2'  Xq 

+ 

-V 

| 2 fi  ' 4 ji  \ -v  ^ ) Wi 

-V 

p ©isj 
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t rq 

V 00 
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CD 
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D 

and  and  are  mode  shapes,  and  i takes  the  values  1 or  2. 


(107) 


2.4.4  Procedure  for  the  Trim  Condition 

The  trim  condition  is  obtained  using  the  iteration  method.  The  pro- 
cedure is  described  as  follows: 

(1)  Advance  ratio  y and  shaft  tilt  angle  are  given  for 
the  wind  tunnel  ejqperiment  condition. 

(2)  Set  an  arbitrary  value  for  the  thrust  coefficient 
CT  as  an  initial  value. 

(3)  Evaluate  the  inflow  ratio  X from  Eq.  106  below  with 

given  y,  a and  C . 

s T 
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(108) 


^ = (U  ta^CXs-  -+ 


C T 


This  is  a fourth-order  polynominal  in  X. 

An  appropriate  root  should  be  chosen  from 
the  resulting  real  and  complex  roots. 

(4)  Using  the  inflow  ratio  A obtained  in  the 
previous  step,  aerodynamic  coefficients 

in  Appendix  A are  evaluated  and  the  equilibri- 
um equation  of  Eq.  103  for  the  no-trim  opera- 
tion or  of  Eq.  105  for  the  moment  trim 
operation  are  solved. 

(5)  Substituting  the  amplitudes  of  8,  and  C , 

lc  lc 

zero  into  Eq.  107  in  the  case  of  out- 

of-trim  operation  (£,  , 8,  and  zero  8. 

lc  Is  lc 

should  be  substituted  for  the  moment  trim 
operation)  yields  a new  thrust  coefficient. 

(6)  Check  the  accuracy  of  the  new  thrust  coef- 
ficient, compared  to  the  old  thrust  coef- 
ficient. If  it  is  not  within  the  accuracy 
specified,  proceed  to  the  next  step. 

(7)  Replace  the  old  thrust  coefficient  with 
the  new  one  and  go  to  step  (3)  for  the 
iteration. 


61 


; 


i 


F 


2.5  Modal  Equations  of  Motion 


In  this  section  the  nonlinear  periodic-coefficient  equations  of  motion, 
Eqs.  34,  35,  36,  and  39  with  aerodynamic  loadings  Eqs.  81,  86,  and  89 
are  linearized  for  small  perturbation  motions  about  the  steady-state 
equilibrium  operating  condition.  Galerkin's  method  is  applied  to  these  equa- 
tions and  modal  equations  are  obtained  in  terms  of  generalized  coordinates. 

The  equations  thus  obtained  are  periodic-coefficient  equations.  To  avoid  dif- 
ficulties in  solving  the  periodic-coefficient  equations,  the  equations  are 
transformed  into  the  nonrotating  coordinate  frame  using  the  Fourier  series  ex- 
pansion, where  most  of  the  periodic  coefficients  are  transformed  into  con- 
stant coefficients.  Stability  calculations  and  frequency  response  analysis 
for  vertical  gusts  are  performed,  based  upon  the  constant-coefficient  equa- 
tions. 

2.5,1  The  Perturbation  Equations  of  Motion 

Bending  motions  in  terms  of  steady  equilibrium  quantities  and  perturba- 
tion motions  in  Eq.  73  are  substituted  into  Eqs.  34,  35,  36,  and  39,  and  the 
equilibrium  equations  are  subtracted.  Discarding  all  nonlinear  products  of 
perturbation  quantities  yields 

w Equation 

'vn.  ur  - + xji)  -4M.XXD.0. Qtn  ^ 

4-  £l(0o  ~Q tw)  Xr'  - 2-wtXlC^  4 {3p4-cTd  ) XT 

- 2vn.n  xrdr 4 [(etf  -+  ei^  uj" 

+ <ev-Eic)6W  v"]"  - tj0  kj'  )'  + 2,  ~ -q  ( ki1  \*v0<X/j  - (-k  rfft)  tr 
+ " ^S^Ca-e.)  - or.J  Gk 


(109a) 


62 


v Equation 
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It  should  be  noted  that  the  single  underlined  terms  are  terms  which 
are  not  included  in  Eq.  40,  used  to  derive  the  natural  frequencies  and  mode 
shapes,  even  though  the  omitted  terms  are  linear.  However,  they  are  retained 
in  the  modal  analysis.  Most  of  these  terms  are  velocity  dependent  terms  re- 
sulting from  Coriolis  forces.  Other  terms,  for  example  (e„,T  4> ) " and 
2 No 

{mx  (r  + r ) d> } ' in  Eq.  109a  where  the  former  expresses  the  centrifugal 
I o 

force  effect,  due  to  the  cross-section  neutral  axis  offset,  and  the  latter  is 
the  centrifugal  force  effect,  due  to  the  cross-section  center-of-gravity 
offset,  are  excluded  in  Eq.  40  because  these  terms  make  the  stiffness  matrix 
nonsymmetrical  as  do  the  Coriolis  force  terms. 

The  double  underlined  terms  are  attributed  to  nonlinear  terms  in  Eqs. 

34,  35,  36,  and  39,  and  include  steady-state  equilibrium  deflection 
effects  producing  structural  couplings.  The  roles  of  these  terms  are  de- 
scribed briefly.  The  term  -mv0®R  in  Eq.  109a  for  the  w-equation  is  due  to  an 
additional  center-of-gravity  offset  resulting  from  the  steady-state  chordwise 

deflection.  The  term  -2mftw'v  is  the  Coriolis  force  effect  from  the  additional 

o 

precone  due  to  steady-state  flapping.  The  rotation  about  the  feathering  axis, 

due  to  the  rigid  pitch,  yields  the  Coriolis  force  -2mv  0 . The  terms 

_ „ OR 

D , R 

-2mf2(w'/  v dr)  and  -2mQ(w"f  v dr)  are  suDplemental  to  the  constant  tension 
r o or 

force  expression  in  Eq.  110,  while  the  original  expression  is  Eq.  37.  The 
2 

term  (-mft  vq0r)  expresses  the  centrifugal  force  effect  due  to  the  center-of- 
gravity  offset  from  the  steady-state  chordwise  deflection.  For  the  rest  of 
the  equations,  short  explanations  will  be  given  if  they  are  different  from 
above. 

In  Eq.  109b  (the  v equation)  the  term  2mf2/r(v'v'  + v'v'  + w'w'  + w'w')dr 

OOOOO 

results  from  the  Coriolis  force  effect  due  to  -2m£2u  in  Eq.  35,  which  includes 
the  geometrical  radial  shortening  effect  from  vertical  and  chordwise  bending 
deflections. 

In  Eq.  109c  (the  4>  equation)  , the  terms  (El  - El  ) [-{w"0  ~ v"}  w" 

F Q O 1W  O 

» ' . ♦ w"}v"]  shows  the  bending-torsion  coupling  identified  by  Mil' 

1 ■■■'.  41).  The  bending  moment  expressions  are  used  to  describe  the 

« . lynamic  loadings  which  cause  torsional  moments  about  the  de- 

. • . < . It  should  be  noted  that  there  is  no  structural  bending- 
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torsion  coupling  in  the  matched-stiffness  rotor  blade  where  El  equals  El  . 

^ c 

For  the  rigid-pitch  equation  of  Eq.  109d,  this  coupling  is  directly  described 
in  terms  of  the  inertia  and  aerodynamic  loadings  and  the  steady-state  blade 
deflections. 


The  steady-state  equilibrium  deflections  and  perturbation  deflections 
may  be  written 


Galerkin's  method  may  be  applied  to  Eq.  109  as  described  in  Subsection 
2.4.2  using  the  natural  frequencies  and  coupled  mode  shapes  in  Eq.  Ill  as 
obtained  in  Subsection  2.2.  Terms  without  underlines  in  the  left-hand  side 
of  Eq.  109  are  converted  into  generalized  masses  and  natural  frequencies,  and 
terms  with  underlines  describe  cross- couplings  between  modes. 


The  modal  equation  for  the  ith  mode,  using  N modes, is  given  by 
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4 Gir  **E  ^<^c/'c*iCf  j + £ { ^«*5  + 

SC  _ ^fcl 

4 h %K  ^ ^ ^ W-</>  c«o</> 

4 ^it*  ?*  4 f*  ^ /^W-  <^1  ] (112) 

where  i = \,  2,  • • ••  Ni  and  5 = lJ~f. 

This  is  the  modal  equation  of  the  individual  blade  where  the  observer  is 
sitting  on  the  blade  and  rotating  about  the  rotor  shaft  with  the  rotating 
angular  velocity  ft.  Therefore,  the  coefficients  of  equations  are  periodic 
in  terms  of  the  blade  azimuth  angle  in  cruising  flight. 

It  should  be  noted  that  Eq.  112  is  applicable  to  the  out-of-trim  con- 
dition in  which  rotor  shaft  tilt  is  prescribed-  However,  no  cyclic  pitch  is 
applied  to  suppress  the  first  harmonic  cyclic  flapping  motion  and  eliminate 
rotor  hub  moments.  Only  the  vertical  gust  is  considered  as  a forcing  function 


» 
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in  the  right-hand  side  of  Eq.  112.  Steady-state  cyclic  pitch  effect  on  the 
moment  trim  condition  and  pitch  control  as  a forcing  function  are  considered 
in  Appendix  C because  they  require  a different  treatment. 


The  various  coefficients  in  the  present  equations  are  shown  in  detail 

in  Appendix  B.  However,  they  will  now  be  briefly  described.  The  generalized 

coordinates  for  the  steady-state  equilibrium  deflections  and  the  perturbation 

deflections  are  expressed  by  q,  and  q.,  respectively.  The  generalized  mass 

k J 

corresponding  to  the  ith  mode  shape  is  1^.  5^  is  the  structural  damping  of 

the  ith  mode  in  terms  of  the  critical  damping  ratio.  Coriolis  force  effects 
are  expressed  by  S„^,  and  is  the  nondimensionalized  rotating  blade 


C..  coefficients  include  single 

ijq 


natural  frequency  at  rotational  speed  SI. 

underlined  terms  in  Eq.  109,  and  e E — ' , E — , E - , E , and 

kljqq  ki  jqq  kljqq  kljqq  kljqq 

E,  . consist  of  double  underlined  terms  in  Eq.  109  which  describe  the 
ki  jqq 

steady-state  equilibrium  deflection  effects  due  to  the  perturbation  blade 


.IS 


■ IS 


motion.  A. .•  and  A. . • are  aerodynamic  damping  terms,  and  A. A?.  , and 
ijq  ijq  ijq  ijq 

are  aerodynamic  stiffness  terms.  Superscripts  IS  and  2C  of  coefficients 


2C 

AT 

ijq 


show  that  these  coefficients  are  coefficients  of  sinSTt  and  cos2f2t,  respec- 

IS-  - - - - is 


tively.  The  quantities  A,  . . — , A,  . ~ • , A,  . ~ ■ • , A,  . A,  , — ,A~~  — , 

kljqq  kljqq  kljqq  kljqq,  kljqq  kljqq 

IS  2S  1C  2C 

A . . — , A , A — , and  A,  . are  aerodynamic  coefficients  related  to 

kljqq  kljqq  kljqq  kljqq  7 2c 

steady-state  equilibrium  deflections.  The  quantities  and  are  aero- 
dynamic coefficients  due  to  the  vertical  gust  where  w is  the  nondimen- 

G 

sionalized  gust  velocity,  w / (f2R)  at  the  center  of  the  rotor.  The  quantities 
■ 2C  IS*  IS 

Gkijq  ' Gkiq  ' Gkiq  ' °kiq  ' Gkiq  ar®  aerodynamic  coefficients  due  to  the 
vertical  gust  and  steady-state  equilibrium  deflections. 


2.5.3  Modal  Equations  in  Nonrotating  Coordinate  Systems 

Equation  112  can  be  solved  in  closed  form  only  when  the  flight  condition 
is  hovering,  in  which  all  periodic  coefficients  vanish.  In  cruising  flight, 
it  is  difficult  to  obtain  closed  form  solutions  because  of  the  periodic  coef- 
ficients. Therefore,  numerical  integration  techniques  are  usually  employed. 
However,  it  is  inconvenient  to  evaluate  the  characteristics  of  such  systems. 
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In  the  present  analysis,  the  harmonic  balance  approximation  method 
(the  nonrotating  coordinate  transformation)  is  chosen  to  evaluate  the  systems 
because  this  method  has  been  shown  to  be  reasonably  accurate  at  advance  ratios 
up  to  0.5  (Ref.  40).  This  range  is  sufficient  for  the  present  purpose,  since 
Eq.  Ill  holds  good  to  an  advance  ratio  of  about  0.4  or  0.5  based  on  the 
assumption  that  the  reverse  flow  is  neglected  in  Subsection  2.4. 

In  the  three-bladed  rotor  which  is  chosen  as  the  rotor  model  in  this 
analysis,  the  perturbation  blade  motion  (t)  of  the  jth  mode  in  the  rotating 
coordinate  system  may  be  written  in  terms  of  nonrotating  coordinates  by  a 
coordinate  transformation  of  the  Fourier  type. 

%j(t)  = %j(t)  4 fj^  CUso(flt)  + ftj  (113) 

For  the  steady-state  equilibrium  deflections,  the  transformation  of 

Eq.  101  is  used.  If  the  jth  mode  corresponds  to  the  flapping  mode,  q?(t), 

1C  is  ^ 

q^  (t) , and  q.  (t)  represent  blade  coning  motion,  longitudinal  flapping 

motion  8,  and  lateral  flapping  motion  8.  , respectively.  It  should  be  noted 
lc  _ _ Is  _ 

that  in  Eq.  101  8,8,  , 8.  , C , C,  , and  Z, , are  constant.  However,  in 
o lc  Is  o lc  Is 

Eq.  113  q° , q^,  and  q^S  are  functions  of  time. 

3 3 3 

Substituting  Eqs.  101  and  113  into  Eq.  112  and  using  a harmonic  balance 
method  results  in  a set  of  ordinary  differential  equations  with  constant 
coefficients  that  can  be  expressed  as 


M(Y}  +[c]{Yj + [K]{Y]  = [p]< 


where  if  N mode  shapes  are  used  {y}  will  be 


^ Cf  ws 


ej°  *t_ic 

OKI  U|>J  OkI  -1 


and  [fa],  [C],  and  [ K ] are  generalized  mass,  generalized  damping  and  generalized 

stiffness  matrices  with  (3N)  by  (3N)  dimensions,  respectively.  [F]  is  a 

generalized  force  row  matrix  with  (3N)  dimensions.  [M] , [C]  force  row  matrix 

with  (3N)  dimensions.  [M ] , [C]  and  [K]  consist  of  submatrices  [m..],  [C..] 

13  13 

and  [k„]»  respectively , which  are  the  3x3  matrices,  and  [F]  of  the  3x1  [f^J 
matrix. 
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The  composition  of  the  matrices  are 


M 


[win] 


— Ouj]--- 


r 


Lf.] 

Cf»] 


owj CVmh]  J 

and  similarly  [c]  and  [K]  matrices  are  composed  as  [m]  matrix 
[m..],  [C..],  [k„]  and  [ f ^ ] may  be  written 
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where  6. . is  the  Kronecker  delta  and  (7  is  /-l. 
il 

6.4  Mathematical  Methods  for  Determining  System  Characteristics 

In  this  subsection,  the  methods  used  to  determine  the  system  character- 
istics are  briefly  described.  These  methods  are  the  eigenvalue  analysis  and 
the  frequency  response  analysis  for  the  vertical  gust  input.  The  eigenvalue 
analysis  is  performed  before  the  frequency  response  analysis  to  confirm  the 
system  stability. 

The  eigenvalue  problem  is  constructed  according  to  Meirovitch  (Ref.  45) 
The  set  of  3N  second-order  ordinary  differential  equations  in  Eq.  114  can  be 
converted  into  an  equivalent  set  of  6N  first-order  ordinary  differential 
equations : 


where 


(x)  = [A]  {*}  + [B] 


(118) 


' } [a]- 

-[Kf'cc]  -UmJ'[k] 

CD  [o] 

CKT'  Up] 

to] 


(119) 


Hence,  the  general  eigenvalue  problem  is  defined  as 

{*}  = Ca]  {*]  (i2°> 

where  fa]  is  a 6Nx6N  nonsymmetrical  matrix.  This  real  general  matrix  eigen- 
value problem  is  solved  by  the  QR  algorithm  described  in  Ref.  46  and  eigen- 
values and  associated  eigenvectors  are  obtained. 

The  frequency  response  problem  is  formulated  next.  Substituting  the 
unit  sinusoidal  vertical  gust  input  w = e1Wt  into  Eq.  114  and  assuming  the 
system  response  with  the  same  frequency  of  the  input  to  obtain  the  particular 
solutions  of  Eq.  114  yield  the  complex  simultaneous  equations 

r h]  ( yj = [p] 

where  [H]  is  a (3N)x(3N)  complex  square  matrix: 


[H]  = -[M](£T+c[c]feHK] 


(122) 


Using  the  Gauss- Jordan  reduction,  the  complex  solution  of  {y}  can  be 
obtained,  and  the  absolute  value  and  the  argument  of  the  complex  number 
gives  the  amplitude  and  phase  angle  of  the  response. 


SECTION  3 


EXPERIMENTAL  PROGRAM 

3.1  General  Description  of  Model 

An  overall  view  of  the  model  rotor  and  mount  is  given  in  Fig.  10.  The 
model  rotor  consists  of  three  blades  with  an  offset  flapping  hinge  at  0.176R 
spanwise  location.  This  flapping  hinge  configuration  was  adopted  to  obtain 
the  desired  flap,  lag  and  torsion  natural  frequencies  with  the  flexure  system 
explained  in  detail  in  a later  subsection.  The  rotor  diameter  is  approxi- 
mately 5 feet  and  the  blade  chord  is  2 inches.  The  blades  have  8 degree 
linear  geometrical  twist.  The  model  rotor  is  an  approximate  Froude-scaled 
dynamic  model  of  a typical  soft-inplane  hingeless  rotor.  The  model  rotor 
solidity  is  0.0704,  which  is  roughly  equal  to  the  solidity  of  most  existing 
full-scale  rotors  in  operation  today.  However,  the  Lock  number  is  2.27, 
which  is  very  low  comparing  to  full  size  rotors.  Additionally,  because  of 
the  orientation  of  the  gust  generator  the  rotor  disc  is  located  vertically 
in  order  to  use  the  lateral  gust  as  the  vertical  gust  in  the  helicopter 
sense.  Blade  characteristics  used  in  the  experiment  are  listed  in  Table  1 
and  Fig.  11.  Detailed  descriptions  of  the  design  considerations  and  the 
resulting  model  may  be  found  in  Ref.  48. 

3.2  Rotor  Blades 

In  the  present  experiment  existing  rotor  blades  are  used,  which  have  a 
NACA  0012  airfoil  section  and  8-degree  linear  spanwise  geometrical  twist. 

To  achieve  the  desired  rotating  natural  frequencies  of  the  blade  lag 
and  torsional  motions  the  flexure  is  placed  between  the  blade  and  hub  as 
shown  in  Fig.  12.  The  offset  flap  hinge  at  0.176R  spanwise  station  produces 
the  equivalent  hub  moment  to  that  of  the  hingeless  rotor.  The  chordwise 
bending  and  torsional  stiffness  of  the  flexure  determines  the  rotating  lag 
and  torsional  natural  frequencies  and  the  spanwise  location  of  the  hinge 
gives  the  appropriate  flapping  frequency.  The  cross-sectional  dimensions 
(width  and  height)  of  the  flexure,  which  was  made  of  17-7  PH  stainless  steel, 
were  chosen  to  obtain  the  desired  stiffness.  In  order  to  evaluate  the 
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rotating  blade  natural  frequencies  from  the  given  stiffness  of  the  flexure 
and  blade,  the  computer  program  based  on  the  analysis  of  the  previous  section 
was  used. 

The  physical  properties  of  the  model  rotor  such  as  mass  distribution, 
stiffness  and  so  forth  were  measured  and  are  shown  in  Table  1 and  Fig.  11. 

Lead-lag  dampers  are  uesed  to  overcome  the  lag  resonance  encountered 
in  the  preliminary  experimental  stage.  The  gravity  force  excited  the  one- 
per-revolution  lag  motion  when  the  rotor  rotational  speed  was  at  the  lag 
natural  frequency  since  the  rotor  was  mounted  vertically,  and  the  structural 
and  aerodynamic  lag  damping  of  the  blade  was  very  small.  It  was  found  that 
the  lag  structural  damping  was  0.2%  of  critical  without  the  lead-lag  damper, 
and  with  the  damper  it  increased  to  7.0%  of  critical. 

3.3  Model  Hub  and  Support 

The  rotor  blade  is  rigidly  attached  to  the  rigid  rotor  hub  with  no  pre- 
cone or  droop.  The  blade  collective  pitch  can  be  changed  at  the  hub  location, 
however,  no  cyclic  pitch  control  is  incorporated. 

The  rotor  hub  is  rigidly  mounted  to  the  rotor  shaft  which  runs  across 
the  tunnel  horizontally  as  shown  in  Fig.  10.  This  is  because  the  rotor  disc 
was  located  vertically  to  obtain  the  vertical  gust  effect  from  the  existing 
gust  generator  at  the  M.I.T.  Wright  Brothers  Wind  Tunnel.  The  rotor  shaft 
is  mounted  3 feet  above  the  tunnel  floor  to  place  the  rotor  in  the  center  of 
the  7 x 10  foot  elliptic  cross-section  of  the  wind  tunnel.  The  rotor  shaft 
is  supported  by  two  trunnions  which  are  connected  to  the  wind  tunnel  balance 
system.  Rotor  shaft  tilting  is  also  possible  to  simulate  the  helicopter 
forward  flight  cruising  configuration.  The  shaft  can  be  tilted  up  to  30 
degrees  forward. 

The  rotor  is  driven  by  a hydraulic  motor  mounted  outside  the  tunnel 
and  a flow  control  valve  in  the  hydraulic  system  is  used  to  vary  the  rota- 
tional speed  between  zero  and  approximately  1200  rpm. 

3.4  Instrumentation 

Each  blade  flexure  was  instrumented  with  strain  gages  in  order  to  detect 
lag  and  torsion  motion.  The  lag  bending  strain  gage  was  located  at  0.128R 
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spanwise  station  (1.1  inch  from  the  cantilevered  flexure  root  which  is  2.4 
inches  long)  and  the  torsion  gage  was  at  0.15R  spanwise  location  (1.7  inches 
from  the  cantilevered  flexure  root) . 

In  order  to  detect  the  flapping  motion  a small  .007-inch-thick  leaf 
spring  was  placed  between  the  blade  and  the  flapping  pin  at  the  hinge  loca- 
tion. The  flapping  hinge  pin,  fixed  to  the  flexure,  is  extended  slightly  and 

slotted  to  provide  a built-in  mount  for  the  root  end  of  the  cantilevered  leaf 

spring.  The  outboard  end  of  the  leaf  spring  has  a sliding  fit  in  a blade- 
mounted  fitting.  The  leaf  spring  provides  a negligible  elastic  restraint  to 
the  blade  flapping  motion.  This  leaf  spring  was  also  instrumented  with  strain 
gages  for  the  measuring  of  blade  flapping  motion  in  terms  of  leaf  spring 
bending  strain.  The  instrumentation  for  flap,  lag  and  torsion  motion  can  be 
seen  in  Fig.  12. 

In  addition  to  the  blade  deflection  signals,  the  rotor  RPM  and  blade 
azimuth  angle  signal  are  available  on  one  of  the  slip  rings.  Once  each  revolu- 
tion the  circuit  on  the  slip  ring  is  closed  and  provides  the  signal  for  the 

i 

RPM  frequency  counter,  while  in  the  recorded  time  history  the  location  where 
the  circuit  is  closed  gives  the  blade  azimuth  angle.  The  rotation  of  the 
shaft  of  the  gust  generator  was  measured  as  the  gust  frequency.  The  gust 
magnitude  was  measured  separately  from  the  rotor  gust  response.  A hot  wire 
was  placed  at  the  rotor  hub  location  and  the  gust  frequencies  were  set  to 
the  same  frequencies  as  those  where  the  rotor  gust  response  was  to  be  measured. 
Thus,  the  gust  magnitude  signal  was  recorded. 

3.5  Rotor  Blade  with  Chordwise  Center-of-Gravity  Shift  Configuration 

The  purpose  of  the  wind  tunnel  tests  is  to  validate  the  theory  developed 
in  the  previous  section  concerning  the  reduction  of  blade  bending  motion  due 
to  a vertical  gust  utilizing  blade  chordwise  mass  unbalance  as  a passive 
control  measure.  One  means  of  accomplishing  the  testing  of  the  blades  with 
different  center-of-gravity  locations  would  be  to  design  and  build  rotor 
blades  which  actually  had  different  center-of-gravity  locations.  This  method, 
however,  would  have  proved  too  costly  and  time  consuming.  Therefore,  it  was 
decided  that  through  the  addition  of  a variable  weight,  placed  at  the  blade 
tip,  the  designed  center-of-gravity  shift  would  be  accomplished. 
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The  scheme  used  to  apply  the  tip  weight  to  the  model  blade  is  shown  in 


Fig.  13.  A stainless  steel  sleeve/tube  is  attached  to  the  blade  tip.  A 
tungsten  slug,  weighing  40  grams,  is  slipped  inside  the  tube.  By  securing  the 
slug  at  any  location  in  the  tube,  an  appropriate  shift  of  the  chordwise  center- 
of-gravity  either  fore  and  aft  was  achieved. 

The  slug  location  allows  for  effective  chordwise  C.G.  variation  (C.G. 
location  averaged  in  the  spanwise  direction)  ranging  from  35%  chord  from  the 
leading  edge  (most  aft  position)  to  18%  chord  (most  forward  position). 

In  order  that  no  more  than  one  blade  of  the  model  be  subjected  to  the 
modification  for  shifting  the  chordwise  C.G.  mentioned  above,  two  dummy  blades 
were  used  in  addition  to  the  actual  blade  having  variable  chordwise  C.G.,  as 
shown  in  Fig.  13. 

3.6  Experimental  Procedure 

Having  finished  the  rotor  static  balance  prior  to  the  wind  tunnel  test, 
the  rotor  blade  track  check  was  performed  with  the  rotor  spinning  at  400  RPM. 
Rotor  blades  were  set  to  0 degree  at  75%  radius.  The  rotor  stability  check 
and  the  structural  integrity  check  of  the  model  rotor  were  also  conducted. 

The  blade  pitch  angle  was  increased  to  4 degrees  and  finally  to  8 degrees. 

The  rotor  rotational  speed  was  then  increased  to  955  RPM  without  tunnel  velo- 
city. The  blade  pitch  angle  8 degrees  and  the  rotor  rotational  speed  955  RPM 
are  the  configuration  for  the  actual  test. 

In  the  next  step  the  rotor  shaft  angle  was  changed  from  0 degree  to  10 
degrees.  The  tunnel  velocity  was  increased  to  30  MPH  and  then  to  60  MPH, 
which  corresponded  to  an  advance  ratio  (y)  of  0.192  and  0.384  respectively. 
Having  again  confirmed  the  rotor  stability,  the  sinusoidal  vertical  gust  was 
applied  to  the  rotor  to  measure  the  frequency  response  of  the  rotor  to  the 
vertical  gust.  The  preliminary  testing  showed  that  the  rotor  response  was 
difficult  to  distinguish  from  noise  at  a gust  frequency  greater  than  500  CPM 
because  of  the  strong  gust  dissipation  at  these  higher  frequencies.  However, 
the  low  frequency  gust  responses  are  of  most  interest.  Thus,  the  wind  tunnel 
test  was  conducted  at  the  gust  frequency  500  CPM  and  below. 

The  sequence  of  the  testing  in  which  data  were  taken  will  be  described 
next.  The  summarized  table  of  test  conditions  is  shown  below. 
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WIND  TUNNEL  EXPERIMENT  CONDITIONS 


TUNNEL 

0 MPH 

30  MPH 

60  MPH 

VELOCITY 

(U=0,  HOVER) 

(p=0. 192) 

(y=0. 384) 

THREE-BLADED  ROTOR 

0,100,200 

0,100,200, 

WITH  STIFF  BLADE 

NO  GUST 

300,400,500 

300,400,500 

TORSION 

AND 

CPM  GUST 

CPM  GUST 

WITHOUT 

TIP  WEIGHT 

FREQUENCY 

FREQUENCY 

APPLIED 

APPLIED 

18%  c 

0,100,150,200 

0,100,150,200 

C.G.  from 

NO  GUST 

250, 300,400, 

250,300,400,500 

L.E. 

500, CPM  GUST 

CPM  GUST 

SINGLE 

FREQUENCY 

FREQUENCY 

"STIFF 

TORSION" 

APPLIED 

APPLIED 

BLADED 

25%  c 

C.G. 

II 

II 

II 

ROTOR 

BLADE 

from  L.E. 

WITH 

35%  c 

C.G.  from 

II 

II 

II 

TIP 

L.E. 

WEIGHT 

18%  c 

C.G.  from 

L.E. 

II 

II 

II 

"SOFT 

TORSION" 

BLADE 

25%  c 

C.G. 

from  L.E. 

tt 

II 

If 

35%  c 

C.G. 

from  L.E. 

II 

• I 

It 

NOTE:  fi  = 955  RPM,  0 =8  deg,  a (SHAFT  ANGLE)  = 10  deg. 

o s 


The  first  test  was  conducted  on  the  three-bladed  rotor  shown  in  Fig.  10. 
The  shaft  forward  tilt  angle  was  10  degrees  and  the  collective  pitch  was  set  to 
8 degrees  at  75%  radius  location.  The  rotor  rotational  speed  was  955  RPM  (100 
rad/sec).  First,  the  hovering  condition  was  achieved  when  the  wind  tunnel  had 
no  velocity.  In  this  configuration  no  vertical  gust  can  be  applied.  The  next 
step  involved  increasing  the  wind  tunnel  velocity  to  30  MPH.  At  this  point 
no  gusts  were  applied  in  order  to  define  the  steady-state  rotor  motion.  Then 
the  vertical  gust  was  applied  to  the  rotor  and  the  responses  of  the  rotor  were 
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measured  at  gust  frequencies  of  100,  200,  300,  400,  and  500  CPM.  Having  com- 
pleted the  30  MPH  sequence  of  the  tests,  the  tunnel  velocity  was  increased  to 
60  MPH  and  the  same  test  sequence  was  followed  as  for  30  MPH. 

The  next  portion  of  the  test  was  with  the  single-bladed  rotor  with  the 
tip  weight  shown  in  Fig.  13.  The  tungsten  slug  was  positioned  in  the  sleeve 
of  the  blade  tip  to  place  the  chordwise  center-of-gravity  at  25%  chord  from 
the  blade  leading  edge.  It  should  be  noted  this  C.G.  location  is  the  averaged 
one  over  the  span.  In  this  case  the  "stiff  torsion"  blade  was  used  (the 
torsional  natural  frequency  is  5.12/rev).  The  rotor  had  8 degrees  of  collec- 
tive pitch,  the  rotational  speed  was  955  RPM  and  the  shaft  tilt  angle  was  10 
degrees.  The  same  test  sequence  was  followed  as  that  of  the  three-bladed 
rotor,  to  measure  the  frequency  response  of  the  rotor  to  vertical  gusts  at 
advance  ratios  0.192  and  0.384. 

The  tungsten  slug  was  moved  forward  to  place  the  chordwise  C.G.  location 
at  18%  chord  from  the  leading  edge,  and  the  gust  frequency  responses  were 
measured.  The  slug  was  then  moved  backwards  to  obtain  the  chordwise  C.G. 
location  35%  chord  from  the  leading  edge,  and  the  frequency  responses  were 
again  measured.  In  the  above  cases  the  "stiff  torsion"  blade  were  used.  The 
rotor  blade  remained  stable  in  any  C.G.  location. 

The  final  part  of  the  test  used  the  "soft  torsion"  blade  (the  torsional 
natural  frequency  is  2.38/rev  with  25%  chord  chordwise  C.G.  location  from  the 
leading  edge)  with  variable  C.G.  location.  The  wind  tunnel  velocity  was  set  to 
30  MPH  and  60  MPH.  The  chordwise  C.G.  location  was  varied  18%,  25%  and  35% 
chord  from  the  leading  edge  at  both  tunnel  velocities,  and  the  vertical  gust 
was  applied  to  the  rotor. 


SECTION  4 


RESULTS  AND  DISCUSSIONS 

4.1  Introduction 

Using  the  theory  derived  in  a previous  section,  the  blade  chordwise 
center-of-gravity  shift  as  a gust  alleviation  method  is  evaluated  and  the 
theory  is  compared  with  the  experimental  results.  Rotor  blade  characteristics 
employed  in  this  analysis  are  based  on  the  wind  tunnel  model  described  in  the 
previous  section  and  detailed  characteristics  are  listed  in  Table  1 and  Fig. 
11.  In  addition  to  the  chordwise  center-of-gravity  shift  as  a principal 
parameter,  the  parameters  examined  in  this  analysis  include  cruising  flight 
at  advance  ratios  of  0.192  and  0.384,  and  blade  torsional  stiffness. 

4.2  Natural  Frequencies  and  Mode  Shapes 

Based  on  the  mass  and  stiffness  distribution  of  the  wind  tunnel  model 
blade  as  shown  in  Fig.  11,  rotating  blade  natural  frequencies  and  associated 
mode  shapes  are  obtained  from  the  theory  developed  in  the  previous  section. 
Mode  shapes  of  the  blade  with  "stiff"  torsional  stiffness  and  C.G.  at  25% 
chord  from  the  leading  edge  are  shown  in  Fig.  14.  Since  existing  stiff  blades 
were  used  in  the  experiment,  bending  in  the  flexure  at  the  root  determined 
the  appropriate  natural  frequencies. 

Non-rotating  natural  frequencies  for  lag  and  torsion  were  determined 
experimentally  with  cantilevered  hub  restraint  and  compare  ' to  the  calculated 
values.  These  are  listed  below.  The  non- rotating  flap  frequency  was  zero 
due  to  the  blade  flap  hinge.  Rotating  blade  natural  frequencies  were  not 
obtained  in  the  experiment. 
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WIND  TUNNEL  MODEL  BLADE 
NON-ROTATING  NATURAL  FREQUENCIES  (Hz) 


LAG  FREQUENCY 

TORSIONAL  FREQUENCY 

EXPERIMENT 

THEORY 

EXPERIMENT 

THEORY 

THREE-BLADED 

9.1 

9.15 

109 

113 

STIFF 

SINGLE- 

BLADED 

TORSION 

7.5 

7.7 

83 

81.9 

WITH  25%  c 

C.G.  FROM 

L.E. 

SOFT 

TORSION 

6.6 

6.9 

38 

39.4 

4.3  Trim  Conditions  and  Steady-State  Deflections 

The  no-trim  steady-state  condition  for  the  three-bladed  wind  tunnel 
model  rotor  described  in  Table  1 is  shown  in  Fig.  15.  Experimental  steady- 
state  deflections  for  flapping  and  lagging  are  also  shown. 

The  theoretical  thrust  coefficient  has  a maximum  point  around  y = 0.1 
and  then  decreases.  This  is  due  to  the  inflow  ratio  decrease  around  y = 0.1, 
which  results  in  an  effective  angle  of  attack  increase.  Large  collective 
lagging  deflection  was  produced  by  the  centrifugal  forces,  on  the  lead-lag 
damper.  The  damper  was  not  scaled  in  the  size  and  weight.  Therefore, 
relatively  large  centrifugal  forces  were  exerted  to  the  damper  which  was 
located  off  the  blade  elastic  axis.  Thus,  the  large  collective  lagging 
deflection  does  not  simulate  the  actual  collective  lagging  deflection. 
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4.4  Frequency  Response  to  the  Vertical  Gust 
4,4.1  Experimental  Data  Reduction 

The  output  of  the  strain  gage  signals  for  the  blade  flap,  lag  and 
torsional  motions  were  recorded  on  magnetic  tape  using  a Hewlett.  Packard  4 
channel  model  3960  tape  recorder.  These  recorded  data  were  fed  to  an  analog- 
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to-digital  converter  (Digital  Equipment  Corp.  computer  PDP-11/40)  and  512 
digital  data  points  were  sampled  with  a 6 msec  sampling  rate.  This  corres- 
ponds to  10.5  data  points  per  rotor  revolution.  After  sampling,  data  were 
digitally  analyzed  by  Fast  Fourier  Transform  (Ref.  49)  with  the  same  machine 
to  obtain  the  frequency  spectrum  of  the  signal. 

Sample  analog  data  from  the  flap  and  torsion  strain  gages  are  shown  in 
Fig.  17.  The  configuration  is  the  three-bladed  rotor  with  tunnel  velocity 
30  MPH  and  vertical  gust  frequency  200  CPM.  From  this  analog  data  the  fre- 
quency spectrum  shown  in  Fig.  18  is  obtained  by  the  method  mentioned  above. 

The  blade  motion  response  to  vertical  gusts  can  be  expressed  in  terms  of  the 
three  frequency  modulated  motions  to,  (ft-<o)  and  (ft+co) , where  ft  is  the  rotational 
speed  of  the  rotor  and  to  is  the  gust  frequency.  With  reference  to  the  blade 
flap  frequency  spectrum  of  Fig.  18,  the  first  peak  from  the  left  is  the 
response  to  the  gust  frequency  to,  called  the  collective  flapping  response. 

The  peak  at  the  nondimensional  frequency  ^ = 0.791  corresponds  to  the  (ft-u>) 
cyclic  flapping  response.  The  big  peak  at  ^ = 1.0  has  ft  frequency  and  comes 
from  the  steady  state  blade  flapping  motion  due  to  forward  flight.  From  the 
present  linear  theory  1ft  steady-state  cyclic  flapping  motion  is  not  affected 
by  the  vertical  gust,  and  this  is  confirmed  by  the  experimental  results. 

The  peak  to  the  right  of  the  ft  frequency  response  is  the  (ft+co)  cyclic 
flapping  response.  The  low  frequency  gust  response  is  of  interest  in 
this  study  and  the  amplitude  of  higher  frequency  responses  such  as  (2ft-co)  and 
(2ft+U))  are  relatively  small.  Therefore,  the  gust  responses  up  to  (ft+co)  fre- 
quency are  considered. 

The  experimental  results  are  compared  to  the  theory  in  terms  of  the 
peak  values  of  the  blade  motion  frequency  spectrum  at  frequencies  co,  (ft-io) 
and  (ft+co) . For  important  responses  such  as  blade  flapping  motion  the  fre- 
quency spectrum  analysis  was  conducted  several  times  on  different  portions 
of  the  magnetic  tape  on  which  analog  signals  from  strain  gages  had  been 
recorded  for  three  minutes.  This  is  a kind  of  averaging  process  because  the 
total  tine  required  for  sampling  is  3.07  seconds  out  of  3 minutes  data. 

During  the  experiment  the  gust  magnitude  and  the  rotor  rotation  speed  are 
fluctuating.  Therefore,  the  rotor  (ft+co)  response  amplitudes  showed  large 
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variation  when  the  summed  frequency  of  the  gust  frequency  and  the  rotor 
rotational  speed  is  close  to  the  blade  flap  natural  frequency,  as  shown  in 
Fig.  20. 

It  is  obvious  from  the  theoretical  analysis  in  Section  2 how  to  obtain 
the  three  frequency-modulated  responses  u),  (ft-u))  and  (ft+w)  . However,  it 

would  be  helpful  to  explain  the  method  briefly  here.  Using  the  rigid  mode 
for  flap  motion,  the  simplified  blade  flapping  equation  with  vertical  gust 
excitation  may  be  written: 

@ 4 P-  -f-  P 4 l C Tl  ) 4 1200  (AO 

+ p.1  ^.ca£J|  p = +1  xoS^dx 

(123) 

where 

x = nondimensional  radial  coordinate 
y = Lock  number 
y = Advance  ratio 

u)0  = Rotating  flapping  natural  frequency 

p 

ft  = Rotor  rotational  speed 
t = Time 

w^  = Vertical  gust  velocity  nondimensionalized  by  blade  tip  speed  (ftR) 
at  the  radial  coordinate  x 

and  if  the  vertical  gust  w is  the  one-dimensional  sinusoidal  gust  in  the 

G 

flight  direction  it  is  expressed  as 

Ulfij  =•  — Gsi)  ~pi  UJ*  csnCQ.'t)  (124) 

The  first  term  in  Eq.  124  represents  the  uniform  component  of  the  gust 
velocity  over  the  rotor  disc,  and  the  second  term  represents  the  gust 
velocity  gradient  due  to  nonuniformity  over  the  rotor  disc.  The  symbol  w 

VJ 

expresses  the  nondimensional  gust  amplitude  at  the  rotor  shaft  location. 
Substituting  Eq.  124  into  the  right-hand  side  of  Eq.  123  yields: 
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[Left-hand  side  of  Eq.  123] 


~ \ /U^-uSt  + 4 { P~4  •+  ^ c^3CCl+uj)t 

+ * ( ~ lKs)|  ^ c^{n~cJ)t  + H.gher  order  Tems  (125) 

Therefore,  due  to  the  vertical  gust  with  frequency  Od,  the  flapping  motions 
0),  (f2-u>)  and  (Q+u))  will  be  excited. 

It  should  be  noted  that  all  gust  responses  in  Pigs.  20  through  23  are 
not  normalized  by  the  gust  magnitude.  The  gust  magnitude  can  be  found  in 
Fig.  19. 

4.4.2  Three-Bladed  Rotor  Gust  Response 

The  three-bladed  rotor  case  provided  basic  information  regarding  the 
overall  response  of  the  rotor  model  to  a vertical  gust  and  also  validated 
that  there  were  no  major  differences  between  the  three-bladed  rotor  and  the 
single-bladed  rotor  in  terms  of  the  vertical  gust  frequency  response.  The 
three-bladed  rotor  gust  response  of  the  experiment  is  shown  in  Fig.  20  for 
the  case  of  30  MPH  (p  = 0.192)  and  60  MPH  (p  = 0.384). 

The  (Q-Hjj)  cyclic  flapping  response  has  a peak  at  a nondimensionalized 
frequency  of  (^)  = 0.2  in  Fig.  20  (a)  and  (b)  since  the  exciting  frequency 
(fHod)  becomes  1.2/rev  and  the  frequency  (ft+Gd)  becomes  close  to  the  blade 
flapping  natural  frequency,  which  is  1.13/rev  in  a vacuum,  and  becomes 
1.2/rev  in  air  due  to  the  aerodynamic  stiffness  caused  by  pitch-flap  coupling 
due  to  steady-state  flap  and  lag  deflection. 

The  theoretical  predictions  of  the  od  collective  flapping  response  at 
both  30  MPH  (p  = 0.192)  and  60  MPH  (p  = 0.384)  are  higher  than  the  experi- 
mental results.  The  collective  motions  are  excited  by  the  uniform  component 
of  the  vertical  gust  and  in  the  theoretical  prediction  this  component  is 
assumed  to  be  the  same  over  the  rotor  disc  as  the  magnitude  at  the  center  of 
the  rotor.  However,  in  the  wind  tunnel  the  vertical  sinusoidal  gust  is  not 
constant  in  the  Y-direction  shown  in  Fig.  9.  At  a distance  of  2 feet  from 
the  center  in  the  Y-direction  the  gust  magnitude  falls  off  to  50%  of  that  of 
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at  the  center.  This  gust  profile  reduces  the  blade  oj  flapping  response. 
However,  in  this  analysis  further  mathematical  modeling  of  the  wind  tunnel 
vertical  gust  profile  has  not  been  conducted.  This  effect  is  not  significant 
for  the  single-bladed  rotor  blade  with  tip  weight  because  of  the  low  Lock 
number. 

The  (fi-u))  cyclic  lag  response  in  Fig.  20(a)  and  (b)  are  higher  around 
(^)  =0.2  because  the  exciting  frequency  (Q— to)  is  0.8/rev  and  close  to  the 
lag  natural  frequency. 

The  (Q-Ha)  cyclic  torsion  motion  is  coupled  with  the  (ft+w)  cyclic  flap 
motion  through  pitch- flap  coupling  due  to  the  steady-state  flap  and  lag 
deflection. 

The  effect  of  the  inclusion  of  the  gust  velocity  gradient  due  to  the 
gust  nonuniformity  is  shown  in  Fig.  21.  Flap  response  levels  are  compared 
with  and  without  the  gust  velocity  gradient.  It  should  be  noted  that  the 
gust  velocity  gradient  has  a significant  effect  on  the  cyclic  flapping 
motion  even  if  the  ratio  of  the  wave  length  to  the  rotor  radius  is  11.1  in 
which  the  advance  ratio  is  0.192  and  the  gust  frequency  u/Q.  is  0.1. 

4.4.3  Single-Bladed  Rotor  Gust  Response  with  Chordwise  Center-of- 
Gravity  Shift 

While  the  initial  wind  tunnel  tests  dealt  with  the  gust  response  of  the 
three-bladed  rotor,  the  most  important  portion  of  the  experiment  involved 
the  single-bladed  rotor  configuration  with  chordwise  center-of-gravity  shift. 
An  objective  of  the  test  program  was  to  evaluate  the  effectiveness  of  tor- 
sional stiffness  variation  in  conjunction  with  chordwise  center-of-gravity 
shifts  in  alleviating  the  gust  response  of  a hingeless  rotor.  The  compari- 
son of  the  theoretical  rotor  responses  with  the  experimental  responses  from 
the  wind  tunnel  test  will  be  discussed. 

4. 4.  3.1  30  MPH  (yi  = 0.192)  Configuration 

The  gust  responses  for  flap  and  torsion  motions  of  the  "stiff"  and 
"soft"  torsion  blade  in  30  MPH  are  shown  in  Fig.  22. 

The  flap  motion  response  with  "stiff  torsion"  blade  are  shown  in 
Fig.  22  (a).  The  (&+ ui)  cyclic  flapping  response  has  a peak  at  the  frequency 
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0.15/rev.  This  is  because  the  (fi-Kd)  cyclic  flap  motion  is  close  to  the  blade 
rotating  natural  frequency  as  explained  in  Subsection  4.4.2.  as  expected, 

the  flap  response  differences  resulting  from  the  chordwise  center-of-gravity 
shift  are  very  small,  because  the  torsional  motion  produced  by  the  C.G.  shift 
was  small  due  to  the  blade  stiff  torsional  characteristics.  In  addition  to 
the  blade  stiff  trosional  characteristics,  the  Lock  number  for  this  single- 
bladed  rotor  is  very  low  at  a value  of  0.951.  This  is  roughly  one-tenth 
that  of  a full-scale  helicopter  and  therefore  the  aerodynamic  forces  on  this 
blade  are  one-tenth  as  effective  in  comparison  with  the  aerodynamic  forces 
on  the  full-scale  helicopter.  Therefore,  the  small  angle  of  attack  change 
due  to  the  torsional  motion  did  not  produce  effective  aerodynamic  forces  to 
suppress  or  excite  the  flap  motion. 

The  torsional  motion  corresponding  to  the  flap  response  in  Fig.  22  (a) 
is  shown  in  Fig.  22  (b) . Since  only  small  variations  are  obtained  in  the 
theoretical  predictions,  the  theoretical  results  of  only  the  25%  chord  center- 
of-gravity  location  from  the  leading  edge  are  shown  in  Fig.  22  (a)  and  22  (b) . 

When  the  "soft  torsion"  blade  is  utilized,  larger  flap  response  level 
differences  due  to  chordwise  C.G.  shift  at  30  MPH  can  be  observed  in  Fig.  22 
(c)  . For  the  (ft+cu)  cyclic  and  u)  collective  flap  response  it  is  clear  that 
the  experimental  response  of  18%  chord  C.G.  location  is  lowest,  that  of  25% 
is  in  the  middle  and  that  of  35%  is  the  highest.  This  tendency  is  the  same 
as  that  of  the  theoretical  prediction  shown  in  Fig.  22  (c) . However,  pre- 
dicted response  levels  for  a given  C.G.  location  are  higher  than  the  experi- 
mental results  of  corresponding  C.G.  location  in  the  (£2+w)  cyclic  and  u) 
collective  flapping  responses.  It  is  difficult  to  describe  the  vibration 
response  reduction  quantatively  from  these  experimental  results.  Correspond- 
ing to  the  flap  response  with  "soft  torsion"  blade  at  30  MPH,  are  the  torsional 
responses  shown  in  Fig.  22  (d) , and  good  agreement  is  obtained  between  the 
theoretical  prediction  and  experimental  observation. 

Lag  motion  responses  of  30  MPH  configuration  are  shown  in  Fig.  22  (e) 
and  (f);  they  are  very  small. 
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4.4. 3.2  60  MPH  (u  = 0.384)  Configuration 

The  gust  responses  at  60  MPH  are  shown  in  Fig.  23.  As  with  the  30  MPH 
configuration,  chordwise  center-of-gravity  shifts  for  the  "stiff  torsion" 
blade  did  not  produce  distinctive  flapping  response  reductions.  Therefore, 
the  theoretically  predicted  25%  chordwise  center-of-gravity  responses  for 
the  "stiff"  blade  are  plotted  to  use  as  a comparison  with  the  experimental 
results.  Overall  agreement  between  theoretical  predictions  and  experimental 
results  is  good.  The  torsional  responses  of  the  "stiff  torsion"  blade  of 
60  MPH  is  shown  in  Fig.  23  (b) . 

The  flap  responses  with  "soft  torsion"  blade  in  60  MPH  are  shown  in 
Fig.  23  (c) . While  there  is  good  agreement  regarding  the  tendencies  for 
the  (f2+u>)  , to,  (fi-co)  flapping  responses,  the  experimental  magnitudes  are 
slightly  higher  than  those  predicted,  particularly  the  (ft+to)  and  to  flap 
responses  of  25  and  35%  chordwise  center-of-gravity  cases.  For  the  ( Q+to) 
cyclic  flap  response  at  the  frequency  oo/Q  = 0.105  the  C.G.  shift  from  35% 


from  the  leading  edge  to  25%  chord  did  increase  the  response  level  slightly 
in  the  experiment  and  no  reduction  can  be  seen  in  the  theoretical  prediction. 
At  the  same  frequency  the  (fi+to)  cyclic  flap  response  is  reduced  by  25% 
experimentally  after  shifting  the  C.G.  forward  from  25%  c to  18%  c chordwise 
station,  when  the  response  level  of  25%  c C.G.  location  is  taken  as  100%. 

In  the  theoretical  prediction  this  reduction  is  20%.  From  the  theory, 
response  reduction  of  the  (fi-Ho)  cyclic  flap  response  can  be  expected  below 
the  frequency  u)/f2  = 0.105  shown  in  Fig.  23  (c)  . This  is  an  important  point. 

A von  Karman  gust  power  spectral  density  in  Ref.  26  and  35  is  given  by: 


(PSD) 


Mi, 
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(126) 


where  a)  is  the  frequency;  L the  gust  characteristic  length;  O the  rms  value 

w 

of  the  gust  velocity;  and  V is  the  aircraft  speed.  Here  L is  approximately 
5000  ft  for  high  altitude  and  from  400  to  500  ft  for  low  altitude  of  the 
full-scale  aircraft.  This  power  spectral  density  scaled  down  to  the  present 
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wind  tunnel  rotor  model  is  shown  in  Fig.  24  where  the  L used  is  40  ft,  V is 
60  MPH  and  the  model  scaling  factor  is  one-tenth.  As  is  shown  in  Fig.  24 
the  gust  power  spectrum  has  a peak  at  the  frequency  ^ = 0.01  and  above 

CO 

q = 0.01  it  decays  very  rapidly.  Therefore,  the  frequency  response  below  a 
frequency  co/ft  of  0.1  is  important. 

Above  the  frequency  ^ = 0.157  the  (ft+ui)  cyclic  flap  response  reduction 
can  be  observed  clearly  in  shifting  from  the  35%  C.G.  to  25%  C.G.,  however, 
there  is  slight  reduction  from  25%  c C.G.  to  18%  c C.G.  in  the  experimental 
results.  Experimental  reductions  of  the  u)  collective  flapping  response  is 
shown  in  Fig.  23  (c). 

Torsional  responses  are  shown  in  Fig.  23  (d)  and  good  agreement  between 
the  theoretical  predictions  and  the  experimental  results  is  obtained. 

Lag  responses  are  shown  in  Fig.  23  (e)  and  (f).  The  theory  predicts 
a lag  resonance  at  the  lag  natural  frequency  and  in  the  experiment  it  was  not 
seen.  It  is  not  clear  why  this  discrepancy  exists.  It,  however,  seems  this 
is  a limitation  due  to  experimental  accuracy  because  the  lag  response  level 
is  very  low. 

4.4.4  Summary  of  the  Vertical  Gust  Response  Analysis 

When  the  rotor  gust  response  is  considered,  the  gust  velocity  gradient 
due  to  the  gust  nonuniformity  over  the  rotor  disc  has  an  important  role  in 
the  blade  cyclic  motion  responses  unless  the  gust  frequency  is  extremely  low 
and  the  gust  uniformity  assumption  is  reasonable. 

The  frequency  response  of  blade  motion  to  a vertical  gust  with  frequency 
co  is  described  in  terms  of  three  frequency-modulated  motions  at  to,  (ft-ui)  and 
(ft+io)  , where  ft  is  the  rotational  speed.  The  responses  (ft— oj)  and  (ft+io)  are 
related  to  the  blade  cyclic  motions.  Therefore,  even  if  the  gust  frequency 
(0  is  low,  the  response  close  to  1/rev  will  be  excited  as  (ft-Ko)  and  (ft-to, 
responses. 

As  is  expected,  the  "soft  torsion"  blade  is  more  effective  in  conjunc- 
tion with  the  forward  chordwise  center-of-gravity  shift  in  alleviating  the 
gust  response  of  a hingeless  rotor  than  the  "stiff  torsion"  blade. 

With  the  "soft  torsion"  blade  the  (ft-Hd)  flap  response  is  reduced  by 
30%  in  30  MPH  and  25%  in  60  MPH  at  the  frequency  ^ = 0.105  when  the  chordwise 
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center-of-gravity  is  shifted  from  25%  chord  from  the  leading  edge  to  18% 
chord  from  the  leading  edge.  This  reduction  is  based  on  the  experimental 
results,  where  the  response  level  of  25%  c C.G.  location  is  regarded  as  100%. 
Theoretical  reduction  is  much  smaller  at  the  frequency  ^ = 0.105  than  that  of 
the  experiments,  but  the  theoretical  prediction  also  show  the  tendency  to 
reduce  the  (£2+w)  cyclic  flapping  below  the  frequency  ^ = 0.1.  Referring  to 
the  von  Karman  gust  power  spectral  density  in  Fig.  24,  which  is  scaled  down 
to  the  present  wind  tunnel  model  and  has  a strong  gust  excitation  below  the 
frequency  ^ = 0.1,  the  forward  chordwise  center-of-gravity  shift  in  conjunc- 
tion with  the  "soft  torsion"  blade  will  be  a good  gust  alleviation  system 
for  the  blade  (Q+oj)  cyclic  flap  vibration.  Response  reduction  of  the  u) 
collective  flap  motion  due  to  the  chordwise  C.G.  shift  was  observed  experi- 
mentally except  at  the  frequency  ^ = 0.105. 

It  should  be  noted  that  greater  effectiveness  of  the  chordwise  center- 
of-gravity  shift  would  be  expected  in  the  full  scale  rotor  since  the  present 
wind  tunnel  model  rotor  has  very  low  Lock  number  (see  Appendix  D) . 

It  should  be  mentioned  that  rotor  rotational  speed  fluctuations  were 
not  considered  in  the  analysis.  However,  these  were  found  experimentally 
to  be  less  than  plus  or  minus  three  percent. 
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SECTION  5 


CONCLUSIONS  AND  COMMENTS 

5.1  Conclusions 

This  study  has  been  devoted  to  the  development  of  a theoretical  model  of 
hingeless  rotor  response  to  a vertical  gust,  and  to  the  theoretical  and  experi- 
mental evaluation  of  the  effect  of  shift  of  blade  chordwise  center-of-gravity 
as  a gust  alleviation  method,  particularly  to  reduce  the  blade  flap  motion. 

Based  upon  the  theoretical  and  experimental  results  in  this  study,  the 
following  conclusions  may  be  stated. 

(a)  Chordwise  center-of-gravity  shift  is  an  effective  and  simple  method 
to  alleviate  hingeless  rotor  gust  response.  In  the  wind  tunnel 
model  blade  with  Lock  number  0.951  and  gust  frequency  0.105/rev 
at  y = 0.394,  it  is  experimentally  confirmed  that  the  (ft+U))  cyclic 
flap  response  is  reduced  by  25%  when  the  chordwise  center-of-gravity 
is  shifted  by  7%  chord  toward  the  leading  edge  from  25%  chord  and 
the  "soft  torsion"  blade  is  employed,  if  the  response  level  of  25% 
chord  center-of-gravity  location  is  regarded  as  100%.  The  theoreti- 
cal reduction  is  20%  at  this  gust  frequency.  Theoretical  prediction 
also  shows  the  tendency  to  reduce  (£2+ol))  cyclic  flapping  below  the 
frequency  c d/ft  =0.1  and  the  response  due  to  atmospheric  turbulence 
will  be  reduced  by  chordwise  center-of-gravity  shift  toward  the 
leading  edge.  Above  the  gust  frequency  aj/fi  = 0.2  at  advance  ratio  of 
0.394  it  is  shown  experimentally  and  theoretically  that  (f2+to)  cyclic 
flapping  is  reduced  by  forward  chordwise  center-of-gravity  shift. 

It  will  therefore  reduce  the  response  due  to  the  discrete  gust 
which  has  higher  frequency  content  than  random  turbulence. 

The  reduction  of  a)  collective  flapping  by  the  chordwise  center- 
of-gravity  shift  is  also  confirmed  experimentally  and  theoretically. 
(f2-to)  cyclic  flapping  response  is  a very  low  level  and  is  not 
significant. 

The  full  size  blade  with  higher  Lock  number  of  order  10  shows 
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even  larger  reductions  of  flapping  response  because  of  greater 
aerodynamic  effectiveness.  (See  Appendix  D.) 

(b)  The  "soft  torsion"  blade  is  more  effective  in  conjunction  with  the 
forward  chordwise  center-of-gravity  shift  in  alleviating  the  gust 
response  than  the  "stiff  torsion"  blade. 

(c)  The  agreement  between  the  experimental  gust  response  of  the  hinge- 
less rotor  and  the  theoretical  predictions  is  good,  particularly 
in  the  important  parameters  of  flapping  and  torsion. 

(d)  The  frequency  response  of  blade  motion  to  a vertical  gust  with  gust 
frequency  w is  described  in  terms  of  three  frequency-modulated 
motions  at  the  frequencies  w,  (ft-co)  and  (ft+oj) , where  ft  is  the  rotor 
rotational  speed.  The  responses  at  the  frequency  (ft-w)  and  (ft+w) 
are  related  to  the  blade  cyclic  motions  and  will  be  excited  even  if 
the  gust  frequency  0)  is  low. 

(e)  The  gust  velocity  gradient  due  to  the  gust  nonuniformity  over  the 
rotor  disc  has  a significant  effect  on  the  blade  cyclic  motion 
responses  unless  the  wave  length  of  the  gust  is  extremely  large 
with  respect  to  the  rotor  diameter. 


5.2  Suggestions  for  Future  Research 

Recommendations  for  the  future  work  before  applying  the  chordwise  center- 
of-gravity  shift  to  the  full  size  rotor  will  be  related  to  two  aspects:  the 
theoretical  analysis  and  the  wind  tunnel  experiments. 

In  the  theoretical  analysis,  the  effectiveness  of  the  chordwise  center- 
of-gravity  to  the  full-scale  rotor  with  Lock  number  10  should  be  studied 


further,  and  it  should  be  confirmed  that  the  center -of-gravity  shift  does  not 
produce  any  adverse  effect. 


In  the  wind  tunnel  tests,  the  response  of  the  rotor  with  center-of- 
gravity  shift  should  be  tested  in  the  very  low  gust  frequency  range  below 
oo/ft  = 0.1.  For  this  purpose  a new  device  to  produce  the  low  frequency  turbu- 
lence and  a more  refined  signal  processing  system  for  the  gust  input  signal 
and  the  response  output  signal  will  be  necessary  to  distinguish  signal  and 
noise. 

The  present  wind  tunnel  rotor  model  might  be  improved  by  the  addition  of 
cyclic  pitch  control.  This  could  be  accomplished  fairly  easily  by  incorporat- 
ing a swash-plate  into  the  model  design,  and  the  necessary  feathering  bearings 
for  this  motion  in  the  rotor  hub. 
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TABLE  1 


DESCRIPTION  OP  THE  ROTOR  BLADE 
USED  IN  THE  WIND  TUNNEL  TEST 


SINGLE-BLADED 
ROTOR  WITH  TIP 
WEIGHT  AND  25% 


THREE-BLADED 

ROTOR 

CHORD  C.G.  LOCATION 

FROM  LEADING  EDGE 

Number  of  blades 

3 

1 

Radius,  R 

2.26  ft 

2.26  ft 

Chord,  C 

2 in 

2 in 

Lock  number,  y (basic  configuration) 

2.27 

0.954 

Solidity,  0 

0.0704 

0.0235 

Collective  pitch,  0q 

8 deg 

8 deg 

Shaft  tilt  angle  in  crusing  flight 

10  deg  forward 

10  deg.  forward 

Lift-curve  slope,  a 

5.7 

5.7 

Drag  coefficient,  CdQ 

0.012 

0.012 

Rotational  speed, 

100  rad/sec 

100  rad/sec 

Built-in  blade  angle  of  twist,  0m, 

TW 

8 deg  (linear) 

6 deg  (linear) 

Elastic  axis 

25%  chord 

25%  chord 

Aerodynamic  center 

25%  chord 

25%  chord 

Precone,  0 

P 

0 deg 

0 deg 

Droop, 

0 deg 

0 deg 

Torque  offset,  eQ 

0 in 

0 in 

Control  linkage  flexibility 

rigid 

rigid 

Basic  Natural  Frequencies 

Lag 

0.823/rev 

0.57/rev 

Flap 

1.13/rev 

1.13/rev 

Torsion 

6. 85/rev 

5.21/rev  (stiff) 

2.38/rev  (soft) 
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UNDEFORMED  COORDINATE  SYSTEMS 
WITH  COLLECTIVE  PITCH  0 


FIG.  5 EQUILIBRIUM  OF  FORCES  AND  MOMENTS 


FIG.  6 DIRECTION  OF  POSITIVE  DISPLACEMENTS  AND  ROTATIONS  AT  NODES  OF 
THE  BEAM  ELEMENT  IN  THE  FINITE  ELEMENT  METHOD 


HUB  PLANE 


(d)  Chordwise  Stiffness 

FIG.  11  STRUCTURAL  CHARACTERISTICS  OF  WIND  TUNNEL  MODEL  BLADE  (STATION 
0.  TO  0.087R  IS  A RIGID  PORTION  OF  THE  ROTOR  AND  FLAP  HINGE 
IS  LOCATED  AT  0.176R) 
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(e)  Torsional  Rigidity 

(1)  Stiff  Torsional  Stiffness 

(2)  Soft  Torsional  Stiffness 


(f)  Angle  of  Twist 


FIG. 


3 SINGLE-BLADED  ROTOR  WITH  TIP  WEIGHT 
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NCNDIMENSIONAL  FREQUENCY  (/rev) 
(a)  FLAP  RESPONSE 


NONDIMENSIONAL  FREQUENCY  (/rev) 

(b)  TORSION  RESPONSE 

FIG.  18  TYPICAL  FREQUENCY  SPECTRUM  OF  THE  VERTICAL  GUST  RESPONSES  OF  FLAP 
AND  TORSION  MOTIONS:  THREE-BIADED  ROTOR  WITH  GUST  FREQUENCY 
u/ft=0.209  AT  p=0.192 
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r 


O--  y=0.192 (30MPH) 


FIG.  19  VERTICAL  GUST  AMPLITUDE  AT  THE  CENTER 
OF  THE  ROTOR  DISC 


2 


NONDIMENS IONAL  GUST  FREQUENCY  u/fi 


(a)  ADVANCE  RATIO  0.192  ( 30  MPH  ) 


FIG.  20  THREE  BEADED  ROTOR  EXPERIMENTAL  VERTICAL  GUST  RESPONSE 

NOTE:  FLAP  AND  LAG  RESPONSES  ARE  EXPRESSED  IN  TERMS  OF  TIP 
DEFLECTION  NONDIMENSIONALIZED  BY  RADIUS.  TORSION  IN  TERMS 


OF  TIP  DEFLECTION. 
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ION  AT  THE  TIP  bOOIMENSIONALIZED 


O I O I O I 


(f2~ (jj)  FLAP  RESPONSE  ^ FLAP  RESPONSE  (Q+oj)  FLAP  RESPONSE 


EXPERIMENT: 


deg  SYMBOLS  AS  ABOVE 


NONDIMENSIONALIZED  GUST  FREQUENCY  m/Q 
(c)  FLAP  RESPONSE  OF  " SOFT  TORSION  " BLADE:  w / « = 2.38  / rev 


FIG.  22  CONTINUED 


EXPERIMENT: 

▲ 

18%  c C.G. 

o 

25%  c C.G. 

0 

35%  c C.G. 

THEORY: 

18%  c C.G 

— 

25%  c C.G 

— 

35%  c C.G 

SYMBOLS  AS  ABOVE 


NONDIMENS IONAL  GUST  FREQUENCY  w/fi 

(a)  FLAP  RESPONSE  OF  "STIFF  TORSION"  BLADE:  u),/ft=5 . 21/rev 

FIG.  23  SINGLE-BLADED  ROTOR  VERTICAL  GUST  RESPONSE  IN  ADVANCE  RATIO 
0.384  WITH  CHORDWISE  CENIER-OF-GRAVITY  SHIFT.  NOTE:  FLAP 
AND  LAG  RESPONSES  ARE  EXPRESSED  IN  TERMS  OF  TIP  DEFLECTION 
NOND IMENS IONAL I ZED  BY  RADIUS.  TORSION  IN  TERMS  OF  TIP  DEFLECTION 
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APPENDIX  A 


EXPRESSIONS  ASSOCIATED  WITH  STEADY- 
STATE  EQUATIONS 

Coefficients  used  in  Eqs.  5.5,  5.8  and  5.10  are  defined  as  follows: 
Inertial  coefficients 

- S.j  cfp+S,,') 2!  - 

■c't  = -(pr  + §p)^  ^ + 

(it) ^ -V  SoG.^  j •vm  (7. - ^ Jvt <Ax 

+ (A.  1) 

Aerodynamic  coefficients 


M 


where 

M 


II 

•0^ 

£ 

+ Mlij 

^Uvv.  Cp 

M«j  - 

ccro^  -+ 

A^-  (></>) 

- 

■ 

■ 

Ms 

e + ri?ce 

/X^~*Sp 

“-j-r 

- Ci- 

(A. 2) 
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I 


4-^Lelj-  +2^Ly-  - gIji  -+(2-  ^ J ^ <5jL ' + 20.  l~&<j 

-e.^] 

^&tj  “ i 1 [ sP^etj  + Ci+  - p - spe„ l'b[j 

-^} 

= T^[2P^^pAy  -pc£j] 

^?;tj  = ^ X pl  Ay  (A>4) 

Mt,  = i X LW?  ^(H|)  AW;*)  + A v«  + Sk-yP*) 

- x1  v>> + -i  p? -w«  _ Go  v«  _ _l  pi v® + >,  e^] 

Ml  = £*!>«!?  - Afxwr'-apav™+Ape„vr] 

Mi.  ■ - M YTfJ  + \ P*  vT  4 - V8(?  - i P* < 

_©„Vcw-xp*e,vrj 

M£>  = i*  [-ifvP-Lfvg  + 4pl0.Vf  ] 

r CK)  r 

-J^^w cTjrfx  , «£’-/'  .«WlV* 
c<;  ^ 

^Qij  ~ I Qs  x<:k)  Wl  VJj  dx  / = f XCK)lfiJ--V}j  cl  I (A.  5) 

cm  ^ ^ 

^ J &S  * / lV  = / Z(K)\fc\J-dj 

q,\  r ^ vl/ 

Wee  - J e,  *<*  WjdX  , W f > - / *•>>  W;  dl 

Vet’  “ / 6*  *wVi.dl  , Vj0')  - / a*’  Vidx 
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A*  = jx«>  W.  Wj  , Cey  -fe-x^Vt- Wj’^ 

Ck)  r 

= J 0s  ^Cc<) 

It  should  be  noted  i and  j take  1 or  2. 


APPENDIX  B 

EXPRESSIONS  ASSOCIATED  WITH  MODAL  EQUATIONS 

Coefficients  used  in  Eq.  6.4  are  defined  as  follows: 

Inertial  Coefficients 

^ = [ [Vw*  + -V  (io+  ^ 

- 2^2*  a^6tw)<5)^  + X (X  XbD~  On  3*  4~~-©tw)Vc©c 

- -a(vn^  ^.©Tw)\rc^: 

- 2.(^vv^.  c^G^'Wc^  4 ^ C^t*~  xSa^1- 

s‘it  ’ -WiVj)rfx 

^ j — "Vn  j /x-UiQt,  Gr^  "VJ i"VJj 

4 2.K  G&.-0rJ  WC^ 

+ H (^)  ^v(6.-  ST»)  vc  wJj  + <*'v  (-^jc^GoC^OrwViV/J^ 
Cfc  = °^(^){cx  c^Qx„  +VcA^-0-rw)C<t-^j  -V^j) 

- «.  C XV/  + xvf) +2(  ^C-t)  {( Wo^. 

+ Vc  'u~Gt“)C-^j  - -*  T„  4’’) 

— . „ J J 

- <feL(Wj  G^^  + Vj-,'^.©TvJ)Jax  <*•» 


where 


*T«  - ( -WvTCdx. 

' X 


Steady-State  Deflection  Coefficients 


4 &j  (yJtMi  -V^0]ax 


*~iti  --  € 


ESlt  i“'[fvc[X(v*V;'+V*A5)di„j 

-•VT^VLVj- V*\7-  \T:]dx 

E^-S,  _ JR?  r 

j l J 'MQc  Cv*  Wj  - \aT^  V-  )d* 

J J (B.3) 

* -g  $ a - [-  V*  Wd  -Wj'  - v*  v;  V/ 

J J 

* vt  f.W/  + v;  vy)*.  - Wi  Wj*  f 1 W(d*„ 

+ M'fW],n 

= S.3  C r r 

*»  3 L 'M  V*  Vi  ©j  - V*  VC®J  - V*  w,  ©^(*4) 

- W^-V' ®;  c^>(2©6)I  _ f _glp  eit  W Ur"-rr"  \ 

J 1 ~ ^ Vj  —c*j 

~ v*-  H'  ***(*&•)  + v*  ^"^ooj  -v  w*V'c<*c2a)j]dr 

Aerodynamic  Coefficients  J 


E^.5,  = R3 


-vtf  - - pO.^? 

- (*«)*$,  + K»b*i 

* A'G’.a.  +w£)b‘£’®;  + (i  (%y  c„,  @L  + ^ C8)V®®| 

+ P & w™ ©-  - 2Ks»-w™ ®j  + n C-£)  cS;®,- ■+  h(|)d&[] 

A^t  -4ft^)(G“’-^-^)] 

A«Jti*iir[^  V^-^%  -2^?H#VC9 

+ + (t!)^  - Ci)vBSUj  - 

■ fS)*«f£®i  + (ll)yAoj  ®'-  4 (■%)  ^®c  + 

- (-&jv,BUj®,  -<£)’  c*®,: 
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- i f.'  (-*■)  C>w  ®; ®J  + xw£ ®c  - ■wj1  + f -wf®; 
-V.f®,  - x^®;  + a'S.^j1’®;  ®j 

4 T H-1  (Pp-VXt>)  Hio-o'  ” ^ ^ ^©0.0 

~ ^ C pp  •+  5t>)  H u/  J 

A'^tt  = i^t  ®C  4 ®L  4 P^J  ®L  - ^-^ei  <©j] 

A'^  +XfAH^'  “P^' 

- ^ %-P  ^o(oj  + ^ - 2Xp.  hS?^  v (^?^SdVVc®j 

-f  ( fr>4  5d)  p L%.  - p (-g ) 4ft,*'  + p C^j  b£ Ij 


-K-&)  4?  fit  + (&)*'&?.<*] 

IS  * ^ 

At)?,  =4v-[2(lW»>(Sj  +2p.L^  -afxC-^B&j 

-ip.(-^r)  B*;®-  - ^-P-ClS)  — ^KCt^)  ^ 

-2yxVe^®,  4A^\  -^pV^Qj  +a^Vi  ®j 


+ 2X|A  L^.  - 2.0„jUl\/dCt)®j  - 20*  p £i<fy  1 

Aa^  = Ar[-p  hSU  -I^F^  +-ip7H^ 


^((ip+So)  V?®j  +i^Cpr^&) 

"ip-  ® A4>:  ^ +(•§ ) fl 1 B4I&  - T (rJK  A ”*  ®c 
+-rP\-%)  Bc2j®i] 

= t V [-  ~ ^Wf  (2)j  - 4;  pX  4 4:  P tfej  ^*i*j 

4 -i  P (-1)  < ®j  4 r h Ck)  ^ 4 4:  K <%) 

-4p  C^-vFd)  H^V  - -ip*  Cp>p4S^p)  Leo^u' 
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r ’ r 


_ ^(^,+40  +i.p.,'e.  . i n'v*"’®, 

/Wjtt  =4rC-^®j  + a^r®j 

+©.^  + a(-%)v^  +*<S)vp£©J  +3(-y„e*!  ®. 

ti (Ss)^" ®c%  - >wW£j  ®4  - -t^®L 

-2Au,vj”’®t  +avV^®t  - +,Vj<*'©:  + &,v5“1®i 

4 ^*(11)  0 0,  ^ 

Asn 4r  l-%if  «t"i  Kx?r«i.  - 

■v  X^®L  + >•  - x y*V®0  + x Vt®;©j 

-V  ^ XJ^  ® ■ - 6W^X)  4 A ®c  <&j  4 4 ^ i 

* i P-vf-A;  -t  ^;Lvcui'  - 

& +^^Sj«t  - k^cV - 

(o')  -Ui)  T— <~°'  -r  T(6W . 

- ®|-  4 ^ 

+ + Nlr^  4 4 ®jl 

~ 2 ^ ^tXW(^4vJi0i  4 ~ *Mr&°‘vj'  " 2P^-4  ^ 

^*■^1  r_2-p-^j‘  ^-2|xV^<SL@j  - 2(J.^^®^0j' 

#.  \ _ 


4}\pTAl£ 

A*.;-  ?<t  = -ivT-l  ilx\3aJ^ 


5 n - K-i  p\*«j  ®c 4 t\$ + i ^“”©,  4 W*}- 


= i r [ wf  1 - vg  + ixvf  - e.  vf  - >C|)  a* 

- (£)  ajf  - *($)am®l  - (jtj  sC3,®c  -Soe**®: 


4 

4- 


A'Sio  © J c^>CX3 


2 * C ^ 


Co) 

el 


4fMM<sx-ap<’] 


* -O-t-M.  C/y 

■ -M-  vK'  +^vt  -c-a^S 

+ (&)^w4  V?®L 
+ c^o<s 

Gi~  “ 

~ ^ P“vjJ^cJ  Cjf»c/y 

-A-W^4-AVec 

-*G.v£  4(fe)A*-  >(^)E4c  4®)ACV 

(13)  'B1 ' *v  A §dG  (5A  +-^ 

4 "z^  cp-^r  - p ■+  Qo ^ ^ <^] 
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- (4  (-|)  - fxSt,  <*=0^ 


Gf  i ~~  2.  & B CSjj'V  "V  |A  (^^>4^1))  \7scJ  /0^w(Vs 

= 4 ~ ^ ®i + ®<i + ^ } &l 

+ 2%  W^®:]  A— ^ 

4 4:^  [ - P-  5 - p <N{X)  c^cy, 

= :k^L~ h^4.  a^njy^. 


G^t 


/2*c 

wL 


_ X 
2. 


i^[-KtoWc  i ^-wj^]  4-WVx 

* [ s^c  + - (-fe)^A<j\:  ~ (-§jwJ^:  - 

~ ^^PCo)®-  + S;  ©^  A-ois 
~ "2  ^ C ~ 6^^  ©C  “*"  ^ ® i-  ~~  "^4  ® l\ 

"iv[-h<  -4-^)a*-4K£)aw 

■v0o^J^(y, 


>!®u 


(B.4) 


y,  «)  r 

AWj  -J'W0^  V;  cU 

Ac<)  r 

A*  -Ja.^v^x. 

"Aj^  */'^t,V*.i;Vj^ 

^L“>  ’ /Am  ^"XVj'cU  ( 
X,,0)  ~ f ' A,.  Vi  AX 


' A*,  -/A.^ajj,/  A% -f A„  f'&Jx 
' A-  =jA"y'*]  Wl  dr^Ai^  , fA„X(“Vt^Vjdl 
, *A?/  - /a,  **  v,  v;  <te/VA«  -/*„  *«H  W/dt 
A^.  - Pv^VjX 
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5©MX«vtax 

e£?  - y^x"'?^,  , wB^j  = jr-^xftJ  W*fe  VjJx  , vp£t  - j5*.xNV*  % 

VB"=  ^ 

e£’  = j ^ , Bw-  fw,  X»\xr  / ^ - s®-  •■*“'% T'Jj  d* 

b2'»^  fl^V^dx 

C^.  = ^ A,  cg’-fc.&QJK  ,c<*>=/chx*>M 

sK)  r 

v^+j  ~ AZ*  ^ egj  oU. , ^cj?  , ^cf05 « Jc,  *00^^ 

PSj  ™ T&l  -=  P^  = [dm^Vl^^I 

dSJ  = Ji^v^x. 

^ = Je^vr^-ax,  p^=  [x^WcXTj^x,  xF^fxCt)V^cVjdx 
v^T  - ^ ^Vs^Wj^  ■*=$*  - f^w,?)ax;  F«J. 

(k.)  r 

- J V]L^  ix,  > fx<">V*  V,  WfdX. 

r-  <X.\  f ' 

= j *co  w^jdix ; 


= { WdVj^X^  <^J  = fx^WV -dl  ,v^j]"f x^tVcy-dx. 
n^J'  -/x^w^VcVj-ax 
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Kiv{Al,  = fe*x^WVj'd< 

uCc'  = { & v,  wj  dx , 5 u'  = j x<Kj  < *J*L  - Vj '<*  X , 

= <J  xCk)  K Vd  tyj-'dx  , * J xCK)  W*  aaTl  Vj'eAx 

r = | x°°  W*V;  Wj  ^ 7 w;WCvUlr'  = J x°°  W*  Vc  ^Jj^x 

LSj  = Vj  ck  , LC*j  =[xwVL  Vjdx,  = /&  *C*° 

hJ-Sj  - f ^ ViUj-dTL  , L%.  = [ Vc$j  <*X  , L?u,u/  - [esX^Vc  V,#dx 

U*ui  " f x°^  Vj'dx  , - [x0^  VL<£j  ax^L^y  =fzCK)W«VL  Vj'AL 

WeL^ - J Xc  W^-(^ ^ dx , ^LS^.  = W«VL ^d\4* no/  = 
w;LCSur  Jx«T<  VLVJ dx  , ^ ^ 
w;L-v.<^  = j"  x°^  W*  \T0  <3jj  di 

Qz)  p 

X7eC  = j ©5  I(t)V,'  d-l  , Vf*  =f  X^  V^dT. , J/^]=  f XCK)  W^Vj dX 
Tr<K}  f r J J 

~Y& i = J Xc°^vjd.x/  vy*-  = 1 0s  ^v v<  Vcd^  Xfcdx 

V*  ~ , J7^=  [OsX^XTe,dX/ vVe^=  |©5XCk)  Y'  Yax 

l^vrdx,  Vyr=[^v;vc^ 


w&(  ~ ( 0S  ^ C<1  y ltff=  | X0^  XKJc dx  , wW-K  = j" x°° ^XFt dx 

'™eJ  = fe.tf^w^cn  .^^je^V^-ax^n 

W*  i :<W  ' Ax  , «WX’  " r©5  Ic"! /XJof  = [ x(WSs T-jV  W: eta. 
5 **  v*  W;*-  , ( xCK)  XjJtVi  dx 
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ft  I fO-1 


J4  = [ w-^^j  ax  J ■x0'1  v^ax  /Xf^t)= 

r=I^^-du 


0<O=  fQ30tc^dX. 


(B.5) 


and 


Am 
Bh 

0,  - 

Dm  = 


- -V  X%S  4l^fi 

<X  a- 


~ Xa 


(X 


(B.  6) 


- (^  4^  7a  + *a) 


T+  -Xa 


In  addition,  the  gust  velocity  gradient  terms  are  expressed  by 
« c 

■ = - ^ , 


oc£t 


= “ -pc  C-^:)  x^Aij 


II 

VV 

yr 

“7*-  Cri)  x 

<=2SC 
Sr  = 

(=r  l/a  = 

<Cc\. 

- 

SC- 

= - TXT  (xx) 

r-3'- 

= "IF Cal) 

rXSc 

~-{r  (&)*•<£% 

- sc 

= " r1  Gs)a-^iii 

(B.7) 


It  should  be  noted  that  the  radial  coordinate  x in  Eq.  B.7  should  be 

• • lc 

integrated  in  the  integrand  of  the  expressions  G^,  and  so  forth. 
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APPENDIX  C 


OUTLINE  OF  EQUATIONS  OF  MOTION  WHEN  CYCLIC 
PITCH  CONTROL  IS  APPLIED 

When  cyclic  pitch  control  is  applied  to  the  rotor  blade,  this  pitch 
control  consists  of  the  trim  cyclic  pitch,  and  perturbation  pitch  control  for 
the  evaluation  as  the  gust  alleviation  system. 

The  trim  pitch  9 is  defined  in  Eq.  5.7,  and  if  it  is  applied  in  moment 
trim  operation,  the  term  0g  in  Eq.  4.18  should  include  this  effect.  Therefore, 

9 in  Eq.  4.19  will  be  rewritten  as 
s 

05  = 0O  + Oic  c^ro^  Q'S  0-tw  (C.l) 

and  the  coefficients  in  Appendix  B should  be  slightly  modified. 

As  for  the  perturbation  pitch  control,  a few  considerations  are  required. 
If  the  rigid  pitch  motion  is  included  and  the  control  linkage  is  flexible, 
rigid  pitch  motion  as  one  degree  of  freedom  is  described  as 

Io  Gl?  Ke  0c  (C.  2) 

where  0 is  the  perturbation  pitch  control.  Therefore,  if  the  control  linkage 
c 

is  flexible,  the  excitation  input  for  the  perturbation  pitch  control  is  only 

the  pitch  control  expression  in  the  rigid  pitch  equation:  for  example,  in 

Eq.  6. Id,  Kfl0  should  be  replaced  by  Ka(0  - 9 ).  Aerodynamic  loadings  due 
w R “ R C 

to  perturbation  pitch  control  are  expressed  by  aerodynamic  loadings  due  to 
rigid  pitch  motion  in  Eqs.  6.1  a,  b,  c,  and  d.  The  final  modal  equations  of 
motion  may  be  written,  modifying  Eq.  6.6 

Cm]  { Yj  +[c]|y}  + [K] (yj  = [pj  jjj.  + [p]  j e*|  (C.3) 

vQisJ 

where 
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APPENDIX  D 


GUST  RESPONSE  OF  THE  FULL-SIZE  ROTOR 


The  preliminary  results  of  the  gust  response  of  the  hingeless  rotor 
with  Lock  number  of  10  are  shown  in  Fig.  D.l.  The  structural  characteristics 
of  the  rotor  blade  are  the  same  as  those  of  the  single-bladed  rotor  with  tip 
mass  used  in  the  wind  tunnel.  The  "soft  torsion"  blade,  = 2.38/rev,  is 

employed.  Blade  density  is  reduced  to  obtain  a typical  full-size  blade  Lock 
number.  The  chordwise  center-of-gravity  is  varied  over  18%,  25%  and  35% 
chord  from  the  blade  leading  edge.  For  the  35%  chord  C.G.  location  the 
classical  flap-torsion  flutter  occurred.  Therefore,  the  results  are  not 
shown  here.  The  gust  magnitude  in  Fig.  19  is  used. 

It  is  shown  in  Fig.  D.l (a)  that  the  center-of-gravity  shift  is  a very 
effective  method  to  reduce  the  flapping  response  due  to  the  greater  aero- 
dynamic forces.  It  is  interesting  that  the  peak  at  0.2/rev  in  the  (Sl-Ho) 
flap  response  disappeared  because  of  the  large  aerodynamic  damping  in  the 
flapping  motion.  The  0)  flap  response  is  higher  than  the  cyclic  flapping 
responses  since  the  gust  effectiveness  has  increased  through  the  Lock  number. 
In  the  (ft-oj)  flap  response  the  response  of  the  25%  chord  C.G.  is  less  than 
that  of  the  18%  chord  C.G.  at  the  gust  frequency  0.4  because  the  phase  of 
the  torsional  motion  is  not  such  as  to  reduce  the  flapping  motion. 
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(a)  Flap  Response 


FIG.  D.l  VERTICAL  GUST  RESPONSE  OF  HINGELESS  ROTOR  WITH  LOCK  NUMBER 
OF  10;  (uym  = 2.32 
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18%  C C.G 


(b)  Torsional  Response 


FIG.  D.l  CONCLUDED 
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APPENDIX  E 


ADEQUACY  OF  HARMONIC  BALANCE  METHOD 

In  this  appendix  the  adequacy  of  the  harmonic  balance  method  used  in 
this  report  is  examined.  The  one-degree-of-freedom  flap  equation  is  employed 
for  this  purpose.  This  equation  is  obtained  from  Eq.  123  as  follows: 

4 It  + (3+  | + ^oao(nt)  4 (3 

= -fii  4 vt-i ^ 

& ' ^ v (E.l) 

~ (£)  coC  otro 
- -£(#)  Kt  c^t;  A^C^rt)  cxro  cat; 

To  compare  with  the  results  of  the  harmonic  balance  method,  the 
numerical  integration  (the  predictor-corrector  method)  is  adopted  to  solve 

O 

the  differential  equation  of  E.l.  The  initial  conditions  are  3=0  and  3=0  at 
t=0.  The  results  show  the  transient  motion  decay  and  the  forced  oscillations 
due  to  the  gust. 

After  applying  the  harmonic  balance  method  to  Eq.  E.l,  the  expressions 
may  be  written: 


O O I 


tit  -2 


-fo-ifO  c&T- 1 


votj  - ^ CCO^t) 


These  equations  are  solved  by  the  same  predictor-corrector  method  with 

o o o 

the  initial  conditions  6 =3,  =3,  =0  and  3 =6,  =3.  =0  at  t=0.  Obtaining  the 

o lc  Is  o lc  Is 

responses  of  3 / 3.  and  3,  at  time  t,  the  flap  response  3 is  calculated 

O 1C  -L  S 

as  follows. 

pet)  = (3o(t)  4 plc(tjc^o (Dt)  4 (E*3) 

It  should  be  noted  that  in  Eq.  E.l  the  periodic  coefficients  higher 
than  one  per  revolution  such  as  sin(2ftt)  and  cos(3fit)  are  included.  However, 
in  Eq.  E.2  they  are  neglected. 

The  results  are  shown  in  Fig.  E.l  to  E. 4.  In  Fig.  E.l  the  gust  response 
with  nondimensional  gust  frequency  0.2/rev  are  shown.  Lock  numbers  were 
chosen  as  10  and  0.954.  The  latter  corresponds  to  the  value  which  was  used 
in  the  present  experiment.  The  flapping  frequency  is  1.127/rev  in  both 
cases.  The  results  from  the  numerical  integration  and  the  harmonic  balance 
method  revealed  very  small  difference. 

In  Fig.  E.2  the  flapping  natural  frequency  is  fictitiously  set  to  3/rev 
to  see  the  adequacy  of  the  expansion  of  the  harmonic  balance  method  written 
in  Eq.  E.3  which  neglects  harmonics  higher  than  1/rev.  The  results  show 
good  agreement  between  the  numerical  integration  and  the  harmonic  balance 
method  when  the  gust  frequency  is  of  the  order  of  0.2/rev. 

The  only  discrepancy  can  be  seen  in  Fig.  E.3  between  the  numerical 
integration  and  the  harmonic  balance  method  when  the  gust  frequency  is 
chosen  as  3/rev  which  is  much  higher  than  normal.  The  3/rev  gust  magnified 
the  motion  in  the  results  for  the  numerical  integration  method  through  the 
3/rev  harmonics  which  are  retained  in  this  method.  In  the  harmonic  balance 
method  these  terms  are  discarded.  Thus,  the  response  of  the  harmonic 
balance  method  is  smaller  than  that  of  the  numerical  integration  method 
for  this  case. 

The  final  figure  E.4  shows  the  results  when  the  flapping  natural 
frequency  of  3/rev  and  the  gust  frequency  of  3/rev  are  chosen.  In  this 
case  the  stiffness  terms  become  dominant  due  to  the  high  natural  frequency 
and  the  aerodynamic  terms  of  3/rev  become  less  effective.  Therefore,  the 
results  show  very  little  difference  between  the  numerical  integration  and 
the  harmonic  balance  method. 
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It  is  concluded  that  the  harmonic  balance  method  used  in  this  study  is 
adequate  for  the  present  rotor  configuration.  However,  there  are  two 
restrictions  when  the  harmonic  balance  method  is  applied  to  general  rotor 
dynamics.  The  first  is  that  the  exciting  gust  frequency  should  not  be  near 
the  same  frequency  as  those  harmonics  which  are  neglected  in  the  method. 

The  second  restriction  is  that  the  eigenvalues  of  the  system  should  be  far 
from  0.5/rev,  1.5/rev  and  so  on,  for  certain  values  of  Lock  number  and 
advance  ratio  (see  Ref.  E.l).  Higher  natural  frequencies  such  as  3/rev  cause 
little  difference  in  the  response  obtained  by  the  numerical  integration  and 
the  harmonic  balance  methods  because  of  the  increased  stiffness. 

Reference. 

E.l  Biggers,  J.C.,  "Some  Approximations  to  the  Flapping  Stability 
of  Helicopter  Rotors",  Journal  of  the  American  Helicopter 
Society,  Vol.  19,  No.  4,  Oct.  1974. 
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APPENDIX  F 

CORRELATION  OF  HUB  SHEARS  WITH  BLADE  FLAPPING  RESPONSE 

In  this  report  the  blade  flapping  response  was  chosen  as  a measure  of 
the  effectiveness  of  the  gust  alleviation  system.  Therefore,  the  appropriate- 
ness of  the  flapping  response  as  a measure  is  examined  by  evaluating  the 
vertical  shear  force  at  the  blade  root,  and  also  the  blade  root  bending 
moment. 

The  vertical  shear  force  at  the  root  is  obtained  by  integrating  the 
aerodynamic  and  inertial  forces  on  the  blade: 

R 

/ ••  \ 

ir 


rR 

S * / (pg  - a)  Jv 


(F.l) 


using  the  force  summation  method,  the  vibratory  shear  force  is  expressed  as 


s 


'■/  \ nHdr 


(F.  2) 


To  evaluate  the  coefficients  of  Eq.  F.2  the  expressions  in  Appendix  B 
are  used,  and  subscript  i is  taken  as  unity.  Mode  shapes  for  i=l  are  speci- 
fied as  W.=W  =1.0  and  V . =V  =$ . =$  =® . =©  =0  for  convenience  in  calculating  the 
i 1 i 1 i 1 i 1 

shear  force.  Subscripts  j are  repeated  from  1 to  4,  which  includes  first 
mode  lagging,  first  mode  flapping,  second  mode  flapping,  and  first  mode 
torsion. 

The  harmonic  balance  method  was  applied  to  Eq.  F.2  to  eliminate  the 
periodic  coefficients,  and  the  shear  force  is  expressed  in  terms  of  co,  (ft-to) 
and  (JHu))  vibratory  shear  forces. 

The  results  shown  in  Fig.  F.l  are  the  calculated  shear  forces  due  to 


obvious  that  the  flapping  may  be  taken  as  a measure  of  the  effectiveness  of  the 
gust  alleviation  system  in  reducing  hub  shears.  In  the  present  wind  tunnel 
model,  the  hub  moment  is  expressed  in  terms  of  the  product  of  the  vertical 
shear  force  and  the  flap  hinge  offset.  Therefore,  the  hub  moment  is  propor- 
tional to  the  shear  force  and  further  evaluation  of  the  hub  moment  is  not 
required. 

The  second  mode  flapping  response  is  shown  in  Fig.  F.2.  It  is 
negligible  compared  to  the  first  mode  flapping  response  when  the  gust 
frequencies  are  as  ”?ow  as  in  the  present  case.  It  should  be  noted  that 
the  shear  force  calculation  in  Fig.  F.l  included  both  the  first  and  second 
mode  flapping  motions. 

As  an  example  of  the  effect  of  higher  harmonic  loading,  the  hub  shears 

due  to  airloads  at  frequencies  up  to  5/rev.  are  presented  in  Fig.  F.3  for  a 

hovering  rotor  of  Lock  number  10.  As  an  expression  of  the  shear  force  the 

2 

nondimensional  form  S/[pA(f2R)  ] is  used  where  S is  the  shear  force,  p the 
air  density,  and  A is  disk  area.  Below  the  gust  frequency  of  1/rev  the 
center-of-gravity  shift  forward  works  very  well  as  a gust  alleviation  system. 
The  peak  around  4/rev  is  due  to  the  second  mode  flapping  response,  and  again 
the  center-of-gravity  shift  reduces  the  response  and  there  seems  no  adverse 
effect  from  the  center-of-gravity  shift. 

It  can  be  concluded  that  the  first  mode  flapping  motion  may  be  chosen 
as  an  indicator  of  the  effectiveness  of  the  gust  alleviation  system. 

It  should  be  noted  that  careful  choice  of  blade  mass  balance  spanwise 
location  and  blade  torsional  frequency  are  necessary  to  avoid  possible 
amplification  of  hub  shears  due  to  higher  harmonic  airloads,  as  discussed 
in  Reference  F.l. 

Reference. 

F.l  Miller,  R.H.  and  Ellis,  C.W.,  "Helicopter  Blade  Vibration  and 
Flutter",  Journal  of  the  American  Helicopter  Society,  Vol.  1, 

No.  3,  July  1956,  pp.  19-38. 
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NONDIMENS IONAL  GUST  FREQUENCY  ^ 

FIG.  F.l  VERTICAL  SHEAR  FORCE  RESPONSE  OF  THE  SINGLE-BLADED  ROTOR  TO 
THE  VERTICAL  GUST  IN  ADVANCE  RATIO  0.384:  "SOFT  TORSION" 

BLADE  Hf)  = 2.38/REV  WITH  CHORDWISE  C.G.  SHIFT 
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FIG.  F . 2 FIRST  AND  SECOND  FLAP  MODE  VERTICAL  GUST  RESPONSE  OF  THE 

"SOFT  TORSION"  BLADE  WITH  18%  c C.G.  IN  ADVANCE  RATIO  0.384 
FLAP  RESPONSES  ARE  EXPRESSED  IN  TERMS  OF  TIP  DEFLECTION 
NONDIMENSIONALIZED  BY  RADIUS 
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